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ACTMO, 

An Agricultural Chemical Transport Model 


By Me H. Frcre,^ C. A. Onstad,^ and H. N. Holtan^ 


SUMMARY 


The use of cheinicals is an important feature of 
modern agriculture, but tire movement of the chemical 
from its point of application constitutes a potential 
hazard* Tire eviiluation of tliis hazard is very difficult 
because the agricultural watershed contains many com- 
plex processes that dynamically respond to weather 
elements. To aid in evaluating tliis hazard, a mathe- 
matical model has been developed and a computer 
program written that describes the movement of chemi- 


cals in and across an agricultural watershed. For each 
storm in series tlie objective is to use tlie model to 
predict tlie concentration of a chemical m the runoff 
water, the total amount carried by tlie runoff water and 
sediment, and tlie location and concentration of the 
chemical remaining on tlie watershed. Tliis report 
describes the rationale used to develop the model as well 
as a couple of tests of tlie model tliat illustrate and 
evaluate its use. 


INTRODUCTION 


Maintaining a higli quality of water resources, con- 
serving land resources, and increasing agricultural ef- 
ficiency are vital ingredients to national sustenance. The 
movement of chemicals and soil from agricultural areas 
can endanger tliese ingredients, 

The Federal Water Pollution Control Act Amend- 
ments of 1972 stress several requirements concerning 
agriculture. These include guidelines for identifying and 
evaluating tlie nature and extent of pollution from 
agriculture lands, and the processes, procedures, and 
methods for controlling it, The Administrator of EPA, in 
cooperation with tlie Secretary of Agriculture and other 
agencies, is charged with developing the necessary 
research to provide the needed information* 


^Soil Scientist, Water Quality Management Laboratory, 
Durant, Okla, 74701. 

^Agricultural Engineer, Soil Conservation Research Center, 
Monis, Minn. 56267. 

^Laboratory Leader, Hydrograph Laboratory, Beltsville, 
Md. 20705, 


Ignorance of essential factors in a complex system 
such as agricultural watersheds can lead to excessive or 
inadequate controls. The result can be economic or 
environmental damage requiring years to correct, The 
factors in a watershed system are interrelated, dynamic, 
and influenced by uncontrollable weather elements. Opti- 
mum management of an agricultural area requires 
predictions of the behavior of the watershed to changes 
in farm management practices. Quantificatioii of die 
potential for any amount of chemical applied to any 
part of a watershed moving to the streams and ground 
water requires a comprehensive sequential accounting in 
time and space for the entire watershed. 

One of the better means for obtaining reasonable 
predictions of complex systems is to perform the 
accounting in terms of a mathematical model. The 
operation of the computational model on modem 
computers permits evaluation of numerous agricultural 
management alternatives and various climatic sequences. 
In addition, testing and development of a model serves 
as a research toot to integrate multidisciplinary projects 
and to reveal gaps in information and understanding. 
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Tile essentials for computing movement of chemicals 
applied to agriculture require a distributed system of 
watershed hydrology, erosion, sedimentation, and chem- 
ical behavior. To establish a continuum for discussion 
and improvement, the USDAHL-74 Model was used for 
watershed hydrology, the Universal Soil Loss Equation 


(USLE) was adapted to compute the erosion, and 
sedimentation on physiographic areas in the watershed. 
A method was developed to compute the interactions of 
agricultural chemicals with the soil -water-plant system 
in the physiographic areas. 


OBJECTIVES 


The objectives for this model are to predict the 
concentration and amount of the chemical in the runoff 
water and on the sediment at the watershed outlet, and 
to predict tlie location and concentration of chemicals 
that are leached and moved spatially through the soil of 
tile waterslied, The model includes such management 
options as the time, the rate, and tlie type of chemical 
applied, changes in crop pattern, and tillage practices. To 
facilitate application, topographic, soil, and land-use 
maps of the watershed are used to estimate model 
parameters. Ultimately, we hope, all watershed informa- 
tion can be obtained from SCS farm plans. We assume 
that published data from the field and laboratory can be 
used fpr estimating tlie chemical interaction with the 
soils and water, 


Considerable research will be required to achieve the 
objectives. This genera! framework constituting initial 
effoits divides tlie model into three submodels: A 
chemical, an erosional, and a hydrological submodcL 
The erosion submodel receives data fiom the hydrology 
submodel while die chemical submodel receives data 
from bo til the hydrological and erosional submodels. 
This separation facilitates tlie interchange of other 
chemical, erosional or hydrological Submodels and modi- 
fications. Objectives of the initial effort were to achieve 
an operating version of ACTMO for certain limited 
conditions. The model must then be tested, modified, 
verified, and extended to major agricultural areas before 
it can be considered a finished product. 


INITIAL MODEL BOUNDARIES 


The watershed is limited to an area tliat can be 
represented by a single rainfall Input (a single gage or a 
composite of several gages), A farm-size area is con- 
venient because it reflects actual management decisions. 
For example, a row crop requiring a large amount or 
number of chemicals can be moved to various parts of 
the watershed. If the scope of consideration is larger 
than a single field, the dilution of the chemical by water 
and sediment from nontreated areas is a major considera- 
tion. The Corn Belt area was used in the initial 
development because it is one of the major areas of 
heavy fertilizer and pesticide use. The initial veision was 
not designed to be applicable to the important agricul- 
tural watersheds located in flat coastal or delta regions 
or containing irrigated or forested land. 

Tlie transport scale is divided into two phases. Hie 
first is chemical application, degradation, redistribution 
by cultivation, and nitrate mineralization and uptake by 
the plants. The second phase is the transportatipn of 
chemicals by water and sediment movement. 

The operation of the model through the first phase 
following application of a nonpersistent chemical is 
adequate to describe tlie potential acute toxicity prob- 


lem. Its operation through numerous phases describes 
the behavior of persistent chemicals that can accumulate 
and pose chrome toxicity or eutrophication problems. 

The initial version is limited to tracing a single 
application of one chemical at a time. The application of 
more than one chemical can be handled by repeating the 
chemical submodel for each chemical. This assumes no 
interaction between chemicals, which is usually the case 
and no direct effect by tlie chemicals on the erosion 
hydrologic behavior of the watershed. FertUizei^ and 
herbicides do affect crop cover and in turn, tlie 
hydrologic and erosion response. This interaction is not 
considered in the initial version. Repeated appUcation of 
me same chemical is important, but programming 
difficulties must be resolved. Salt movement and salinity 
are additional problems not included in this veision. In 
fact, the characteristics of these problems are sufficiently 
different tliat an entirely new chenucal section may be 
required. 

The erosion submodel is limited to rill and interrill 
components of sheet erosion. However, the different 
characteristics of rill armonng and preformed rills are 
not considered. Local deposition is predicted througliout 
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a waterslied, but differences in erodibility and texture of m the watershed. Therefore, chemical considerations 

tJie deposited material aie not considered for subsequent limit a single crop for each zone because average vaiues 

storms, Sediment from areas such as gullies, roads, and for dissimilar crops will not provide adequate infonna- 

construction sites is not accounted for, althougli such tion for chemical movement. Only two layers per zone 

sediment may be significant in downstream chemical are physically defined and provide useable data. Two 

transport. Tlie interactions of fertilizers, soil condi- other layers are empirically defined as storage from 

tioners, and otlier cultural activities witli the soil recession flows. The subsurface movement of water 

properties known to affect erosion are not considered, laterally from each zone is routed across the alluvium or 

nor is tlie time rate-of-change of properties consideied. directly into tiie channel. Snowmelt is computed in the 

The hydrology submodel computes average values for model for all precipitation increments identified as snow 

evapotran- piration, infiltration, and runoff for each zone in the input data. 

THE GENERAL FRAMEWORK 


As previously pointed out, ACTMO is divided into 
three sections to provide a means where alternate 
chemical, erosional, or hydrological submodels could be 
easily substituted. Figure 1 summarizes the general 
framework of the model in a flow diagram while table 1 
describes the information transmitted in the input and 
output files of ACTMO. 

The USDA Hydrograph Lab model was adopted as 
tile hydrologic submodel because it was operational, was 
calculating most of tlie variables considered necessary 
for a transport model, and utilized parameters readily 
available from soil suivey data. The hydrology submodel 
was developed from data on a variety of watersheds 
across tlie country and, therefore, is potentially broad 
based. The eiosional and chemical submodels were 
developed specifically for this model, combining many 
previously published concepts. 



I OUTPUT ) 
\ OF / 
CHEMICAL 

behavior 


Figure l.-Flow diagram for the Input-output files in ACTMO. 


The processes involved m chemical transport indicate 
file following factors are important. The soil moisture 
and temperature influence the breakdown of many 
chemicals and the mineralization of nitrogen, but nitro- 
gen uptake is expected to be correlated with transpira- 
tion. The leaching efficiency of water tends to be lower 
as gravitational water percolates tlirougli the soil. The 
amount of runoff is obviously an important variable as is 
the relative amount of lateral flow for mobile chemicals, 
Therefore, the hydrology submodel supplies the chemi- 
cal submodel the following variables: The temperature, 
and for each layer in a zone the evapotranspiration and 


Table L~Description of the information transmitted in 
the input-output files of ACTMO 

File Information 


Weather . . . 

Ralnfail, temperature, pan evaporation 

Watershed . . 

Zone sizes and cascading, soil hydrologic proper- 
ties, crop pattern and characteristics, tillage 
operations. 

Storm .... 

Length, slope, and cascading of zones, runoff, 
infiltration, and rainfall intensity. 

Soli 

Soil texture. 

Erosion .... 

Compartment areas, enrichment coefficient, soil 
texture, infiltration, runoff, erosion, and deposi- 
tion. 

Hydrological 

Cultivations, soil properties, time-averaged tem- 
perature, soil moisture, and ET, lateral and 
vertical flow, soil water when the storm starts. 

Chemical , . 

Date, amount, and type of chemical application; 
treated area; soil chemical interactions; and 
critical levels, 

Output of 
Chemical 
Behavior. . . 

Chemical lost in runoff and erosion. 
Concentrations in runoff, ' 

Chemical left in each compartment. 

Location and concentration of peak in each 
compartment. 
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soil moisture averages between storms, the infiltration 
before runoff starts, infiltration during runoff, volume 
of runoff, die lateral and vertical flow rates for each 
zone, and tlie soil moisture in each layer of each zone at 
the start of the storm. 

Chemicals adsorbed to tlie soil are subject to erosion 
as a transport mechanism. The chemical concentration in 
the soil varies witli deptli because of the application 
method and the amount of leaclung. Therefore, both rill 
and interrili or surface erosion are calculated. Finally, a 
decrease in the transport capacity of the runoff can 
cause deposition of the sediment, Chemicals are usually 
preferentially adsorbed by the finer sediment fractions 
which are tlie last to be deposited and, therefore, the 


composition of the sediment which is enriched with the 
fine material is calculated. In summary, the erosion 
submodel supplies the chemical submodel with estimates 
of interrili erosion, rill erosion, area and depth of rills, 
amount of deposition, and the enrichment of the eroded 
sediment relative to the clay fractions of the original 
soil. 

The Universal Soil Loss Equation was modified to 
compute the erosion and deposition information for the 
chemical submodel. The rainfall factor in this equation 
was adjusted from an annual average to a storm basis by 
the rainfall intensity, tlie amount of runoff, and the 
peak rate of runoff variables supplied by the hydrology 
submodel. 


THE WATERSHED TRANSFORMATION 


On a multicropped watershed one must anticipate 
til at only part of tlie watershed might be treated with a 
chemical. Some means must be used to transform the 
watershed defined by physical characteristics, such as 
crop fields, into computational units. A zone, the 
computational unit for the hydrology submodel, is 
constructed by grouping fields of tlie same crop, or 


major soil areas, or land-capability classes or whatever 
physical feature is recognized as important. The zones 
typify the elevation sequence of uplands, hillsides, and 
bottomlands. The percent of each zone contributing 
surface runoff to lower zones is estimated by observa- 
tions from the soils and topographic maps. The average 
downslope dimensions of each zone is tlie length of 


ZONES 



Figure 2. “A three-dimensional transform of the waterslied into computational units. 
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overland flow in the zone. The results are illustrated in 
figure 2 where foui zones with different slopes are each 
divided into four compaitments creating seven tubes. 

Each tube is the area between adjacent streamlines 
from flow net tlieory. The flow of water and sediment is 
assumed to occur along the line of steepest descent, 
representing flow over an isotropic surface. This provides 
a situation analogous to the natural sequence of flow 
while maintaining the distributed effect. The size of the 
incremental areas within a zone or compartment is 
calculated from the average zone length and the fraction 
of the zone area that contnbutes directly to the channel. 
Each compartment is intended to be a homogeneous 
unit for computing infiltration, evapotranspi ration, ovep 
land flow, and net erosion or deposition. Soil and water 

SUBMODEL 

Tliis section provides a detailed description of each 
submodel. Such a description is not essential for an 
overview, but a model, like a chain, is no better than its 
weakest link. In tliis section all of tlie lelations and 
assumptions are laid bare for examination. 

The description of each submodel is essentially an 
independent paper. The section on the hydrology 
submodel is a recent revision of “USDAHL-70 Model of 
Watershed Plydrology’* {2G)f^ Tlie erosion-deposition 
submodel was developed by C. A. Onstad while M. H. 
Frere developed the chemical and nitrate submodels. 

The Hydrology Submodel for ACTMO 

The USDAHLr74 Revised Model of Watershed Hy- 
drology^ 

Hydrologic Grouping of Soils and Land Use 

Soils on each watershed are grouped by land capa- 
bility classes, as illustrated in figure 3, to form hy- 
drologic response zones for computing infiltration, 
evapotranspiration, and overland flow. Zones indicated 
for each watershed typify the elevation sequence of 
uplands, hillsides, and bottom lands in these areas. Zones 
are always numbered in a downslopc order because 
computations assume that some of the runoff will 
cascade over successive zones. For example, on the 
Coshocton watershed an estimated 60 percent of Zone 1 
runoff flows directly onto the alluviums of Zone 3, and 
only 40 percent cascades to soils of Zone 2. 

^Italic numbers in parentheses refer to Literature Qted, 
p. 37. 

^USDA“ARS Plant Physiology Institute Report No, 4 by 
H. N, Holtan, G. J. Stilner, Wm. H, Henson, and N. C. Lopez, 
1974. 


lost from one compartment are considered to be 
additions to tlie top of the adjacent downslope compart- 
ment wheie siniOar calculations are made. This process is 
continued down to the watershed mam channel. 

Figure 2 also illustrates tlie vertical dimension that is 
considered. The hydrology submodel treats the first 
layer as the ‘‘A” hoiizon oi plow layer and the second 
layer as the remainder of tlie potential loot zone. The 
other two layei-s are defined empirically iis storage from 
recession outflows. The chemical submodel uses the 
same layers to compute the chemical interactions with 
tlie soil. Althougli coefficients of adsorption, dispersion, 
and breakdown could be handled as continuous func- 
tions of depth, the degree of lefinement is limited by the 
degree of hydrological definition of the soil. 

DESCRIPTION 

The watershed description data published by the 
Agricultural Research Service field stations {48) give 
percentages of area of major soils and of land-capability 
classes. Moisture-tension data (table 2) published by the 
USDA Hydrograph Laboratory (22) or the Soil Survey 
Investigations Reports of the Soil Conservation Service^ 
provide soil profile descriptions for computing weighted 
averages of hydrologic capacities (5) of soils in each 
zone. Specific gravity was estimated at 2.65 for all soils 
in computing total porosities. Moisture at 0.3 bar 
tension was used for field capacity in dividing total 
retention capacities (S) into moisture freely drained by 
gravity (G) and moisture drained by vegetation (AWC^ 
Tlie lower limit of AWC is assumed to be tlie moisture at 
15 bais tension. Table 3, by C. B. England, is a summary 
of a great many samples giving S, G, and AWC m percent 
volume of textural horizons. Percentages are applied to 
depths of **A’' horizons or to plow depths in nonlayered 
soils, and to depths of maximum root penetration as 
miglit be limited by hardpans, gravelly layers, or other 
textures unsuited for root growth. The input format of 
the model accepts G and AWC in percentages by volume 
and applies them to soil deptlis to compute volumes. 

Evapotranspiration 

Evapotranspiration potentials are estimated by coef- 
ficients applied to publislied pan-evaporatipn data. The 
method is a combination of techniques developed by 
Mustonen and McGuiness (.?2), Jensen (25), and 
Pruitt {38) adapted for use with cardinal temperatures 
presented in table 4 in estimating plant growth. Figure 4 
is a plotting of XGI, cardinal temperature functions 

^SoU Conservation Service, Soil Survey Investigations Report 
Series by States, U S, Dept. Agr., Soil Conserv. Serv. 
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Table 2 -Profile description and physical proper ties- Cecil Gritty Sandy Loam 


Location Blacksburg, Virginia, Brunswick County, Rocky Run Branch Watershed, located in pasture 75' SW of 
fence, 400 yds. W of county road #656. 

Vegetation and land use: Pasture. 

Topography. Gently sloping. 

Drainage Well-drained 

Parent material’ Granite gneiss and schist. 

Described and sampled by. K. Fussell, J. Williams, and J. B. Burford. 




Available moisture 

Horizon 

Description 

Per inch of 
soil 

In horizon 

Ap 

0 to 6 inches. Dark yellowish brown (lOYR 4/4) gntty sandy loam, very friable; weak 
medium granular structure; many fine roots, medium acid; few small peds of red (2,5 YR 

4/6) brought up from lower horizons by plowing; abrupt wavy boundary. 

Inches 

o.io' 

Inches 

0.60 


6 to 15 inches. Red (2.5 YR 4/6) clay with few fme mottles of yellowish red (5YR 5/8), 
some ped faces are coated with yeUowish red (5YR 4/8) clay, which is apparently result 
of root and worm action; firm; moderate medium subangular blocky structure, thin 
continuous clay films; few mica flakes; common fme roots; occasional quartz rock less 
than 1/2^^ m diameter; medium acid; clear wavy boundary. 

.07 

.63 

B23t 

IS to 31 inches. Red (2 SYR 4/6) clay loam with common fine faint mottles of yellowish 
red (SYR 5/8) and common fine distuict mottles of yellowish brown (lOYR 5/8); moderate 
medium and fine subangular blocky structure; friable; thm clay film on most ped faces; 
common mica flakes, 10 percent gnt by volume, few fine roots; strongly acid; clear wavy 
boundary. 

.06 

.96 

B31 

31 to 42 inches. Red (2. SYR 4/6) clay loam with many medium faint mottles of yellowish 
red (SYR 5/8) and few medium distinct mottles of strong brown (7. SYR 5/8) and reddish 
yellow (7,5 YR 6/8), friable; moderate to weak medium and fine blocky and subangular 
blocky structure; patchy clay film; many mica flakes; 10 percent grit by volume; strongly 
acid , clear smooth boundary. 

.13 

1.43 

B32 

42 to 56 inches. Mottled colors of red (2. SYR 4/8), yellow red (SYR 5/8), reddish yellow 
(7.5 YR 6/8), very pale brown (lOYR 7/4) highly weathered granite gneiss and schist with 
few pockets of light clay loam, friable, weak fine and medium subangular blocky structure; 
many mica flakes; shmy and greasy to the feel, 10 percent grit by volume; few patchy clay 
films; very strongly acid; gradual wavy boundary. 

.18 

2.52 

C 

56 to 71 inches plus Mottled colors of yellowish red (SYR 5/8), red (2 SYR 4/6), reddish 
yellow (5YR 6/8), and very pale brown (lOYR 7/4) highly weathered granite gneiss and 
schist; silt loam with 10 percent gnt by volume. 

.18 

— 


Horizon 

1 Depth 

Bulk 1 
density 

i 

Moisture as percent of dry weight at respective atmospheres of 
tension 



0 1 

0.3 0.6 

3 15 


Inches G /cc. 

Ap 0-6 1.72 

D22t 6-15 1.33 

B23t 15-31 1.32 

B31 .’ 31-42 1.41 

B32 42- 56 1.29 

C ' 56-71-h 1.29 


Percent 

Percent 

Percent 

Percent 

Percen t 

15.24 

12.35 

10,30 

9.21 

6,34 

26.90 

23.76 

22.22 

21.53 

18.28 

32.02 

29,26 

28.67 

27.27 

24.58 

29.89 

28,02 

25,87 

21.68 

18.62 

35.15 

26.66 

19.52 

16.49 

13.02 

36.36 

24.46 

26,29 

15.78 

10.69 
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Figure 3. -Hydrologic zones of land-capability classes. 


Table J, -Hydrologic capacities of soil texture classes 

[Adapted from “Land Capability: A Hydrologic 
Response Unit in Agricultural Watersheds/* by C, B. 
England, Res. Soil Scientist, USDA Hydrograph 
Laboratory, Beltsville, Md. (USDA ARS 41-172, 
SepL 1970)] 


Texture class 

S* 

G^ 

AWC^ 

AWC/G 


Percent Percent Percent 


Coarse sand 

... 24.4 

17.7 

6.7 

0.38 

Coarse sandy loam .... 

. . . 24.5 

15,8 

8.7 

0.55 

Sand 

. . . 32.3 

19.0 

13.3 

0.70 

Loamy sand 

37.0 

26.9 

10,1 

0.38 

Loamy fine sand 

... 32.6 

27.2 

5.4 

0.20 

Sandy loam " 

. . . 30.9 

18.6 

12.3 

0.66 

Fine sandy loam 

... 36.6 

23.5 

13.1 

0,56 

Very fine sandy loam , 

. , . 32.7 

21.0 

11.7 

0,56 

Loam 

... 30.0 

14.4 

15,6 

1.08 

Silt loam 

... 31.3 

11.4 

19.9 

1.74 

Sandy clay loam 

... 25,3 

13.4 

11.9 

0.89 

Clay loam 

. . . 25.7 

13.0 

12.7 

0.98 

Silty clay loam 

... 23,3 

84 

14.9 

1.77 

Sandy clay 

. . . 19,4 

11.6 

7.8 

0.67 

Silty clay 

, . . . 2L4 

9.1 

12.3 

1.34 

Clay 

. ... 18.8 

7.3 

11.5 

1.58 


^ S = Total porosity - 15 bar moisture percent, 
^ G = Total porosity - 0.3 bar moisture percent 
^ AWC = S minus G, 


(current temperature - lower cardinal temperature)/ 
(optimum cardinal temperature - lower cardinal tempera- 
ture) for corn and alfalfa. Cardinal temperatures are the 
optimum and lower or minimum temperatures for crop 
grow til. Wlien current temperature exceeds the optimum 
cardinal limit, the plant is assumed to suffer and the XQl 
function is set less tlian 1.0 by one-half of the amount of 
exceedence. If current temperature is equal to or less 
tfiaii the lower cardinal limit, the XGI is set equal to a 
fixed minimum of 0.1 , 

The growth index for a given crop is then computed 
from the XGI curve for that crop, as illustrated in figures 
5 and 6 for corn and alfalfa respectively, Tlie crop 
growth index, GI, follows tlie XGI if no cultural 
practices intervene. However, cultural practices that 
reduce tlie foliage also reduce evapotranspiration; hence, 
GI for the crop is reduced following plowing, planting, 
cultivating, or harvesting if vegetation is abolished or 
significantly reduced by the operation. (This excludes 
minimum tillage practices,) After each practice, the crop 
is assumed to recover in a specified number of weeks at 
which time the GI curve rejoins the XGI curve. In the 
case of tumplowing, the soil is inverted and left in an 
aerated condition; hence, the GI value is set at 1.0 and 
assumed to settle back to a bare fallow condition of low 
ET in two weeks wliere it remains a constant at 0.1 until 
planting of a crop. 
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Cardinal temperatures can be obtained from the 
literature for most crops. Table 4 is a compilation of 
upper and lower cardinal temperatures for some of the 
more common crops togetlier with literature citations. 
GenerciUy, tlie local reputation of a crop for seasonal 
growth can be drawn upon to approximate the cardinal 
limits; if a grass is known to go dormant at 


midsummer, the associated temperatures indicate the 
upper limit. Winter wheat would be expected to have a 
lower limit near to freezing. The evapotranspi ration 
potential can then be estimated for any given day as 

ET - GI • k • Ep • [(S-SA)/AWC] "" [1] 


Table 4, -'Cardinal temperatures for growth of some common cropSf degrees Fahrenheit 


Crop — 


Cardinal temperatures 


- Reference source 


Mimmum 

Optimum 

Maximum 


Oats, rye, 
wheat, barley. ^ 

Sorghum, 
melons . . * . * 

32-41 

59-64,4 

77-87.8 

87.8-98.6 

87.8-98.6 

111.2-122 

Parker, *N. W. 1946. Environment factors and their control m 
plant environments. Soil Sci. 62:109-119, 

Corri 

Peas 

Citrus fruit ... 

50 

40 

55 



Chang, J. H. 1968. CUmatc and Agriculture: An Ecological 
Survey. Aldine PubI, Co., pp, 75-76 

Oats. 

40 



Wiggms, S, C. 1936, The effect of seasonal temperature, . . 

Agron. L 48. 21-25. 

Corn. . 

50 

86-95 

113 

Mcycr, B. S. and D, B. Anderson, 1939. pp, 589*591, P. van 
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Figiiru 5, -Relative growth (GI) curve for corn as modified by tillage. 
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( igurc 6 Keldtivc growth (GI) curve for alfalfa as modified by cutting 


wlu're 

I-Tb cvapDtranspiratiun poientul m inches per 
day, 

(jI is growtli index of crop m percent of maturity, 
kis ratio ol (>1 to pan evaporation, usually 
I 0-1.2 for ^hort grasses, 1.2-1 6 for crops up to 
slioulder-hejglit, 1 6-2 0 lor lorcst, 

t-p b pan evaporation in inches per day. 

S b total piuosity. 

SA h available porosity, 

AWT IS fKiri>sity drainable only by cvapotranspira- 
tion, 

and \ IS set ecjual to 0.5 for sands and 1.0 for finer 
textures 

Input to the niodel now includes 52 weekly averages 
ot air temi>eraiure in place of the Gl curves required for 
each crop in Model-70 The use of leinperature is 
designed to indivi dualize plant growth estimates for each 
vear. 


Dr Wadleigh, in mathematical analyses of plant 
growth data (unpublished), demonstrated that free water 
(above field capacity) dampens and inhibits plant 
growth, and subsequently ET, by exclusion of oxygen 
needed for plant growth. As illustrated schematically in 
figure 7, growth and, consequently, ET decreases from 
optimum at field capacity to zero at soil salination. 
Thus ET increases from wilting point to an optimum a t 
field capacity and diminishes to zero at saturation. 
However, free-water evaporation mcreases from zero at 
field capacity to a value equal to pan evaporation nt soil 
saturation. Tins separation of evaporation and transpira- 
tion becomes important in chemical transport studies. 

Infiltration 

Infiltration capacity was expressed by Holtan(7(5) as 
an cxiiaiistion phenomenon convergent upon some 
constant rate* 

f=3 • Sj' '* + fg rj. 
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where 

f is infiltration capacity in inches per hour; 
a is tlie infiltration capacity m inches per hour per 
(inch) of available storage (index of surface- 
connected porosity); 

Sg is available storage in the surface layer (the “A” 
horizon in agricultural soils) in inches water 
equivalent; 

and fc is the constant rate of infiltration after pro- 
longed wetting (associated with capillary flow or 
with an impeding stratum) in inches per hour. 

The G1 of equation [1] is also used as a seasonal 
factor on tlie vegetation parameter ‘*a” in the infiltration 
equation [3] : 

+ fc [3] 

Infiltration and rainfall excess are computed for each 
soil zone by comparing observed rainfall to the infiltra- 
tion capacities computed by equation [3], In solutions 
of equation [3j » the volume of infiltration cannot 
exceed the volume of water available at the soil surface 
(including detentions of overland flow), nor can it 
exceed the product of the time increment and the 
average rate, (fi+f2)*t/2. Infiltration capacity de- 
pends on the amount of empty space with Sa diminished 
by infiltrated water but recovered by drainage or 
evaporation to tlie limit of G, and by ET to the 
additional limit of AWC. 

Values of are sliown In table 5 for various land 
uses in agriculture. Tliese were evaluated at plant 
maturity as the percentage of the ground-surface area 



Figure 7, “Relative growth as modified by free water (matrix 
suction less than .25 bars) and moisture deficit (matrix suc- 
tion greater than .25 bars). (After Wadleigh, unpublished.) 


Table 5. - Tentative estimates of vegetative parameter 
in infiltration equation 


f=Gf* a - 

' a c 


Land use or cover j 

Basal area rating^ 

Poor 

condition 

Good 

condition 

Fallow^ 

0.10 

0.30 

Row crops 

.10 

.20 

Small grains 

.20 

.30 

Hay (legumes) 

.20 

.40 

Hay (sod) 

.40 

.60 

Pasture (Bunch grass) . 

.20 

.40 

Temporary pasture (sod) 

.40 

.60 

Permanent pasture (sod) 

.80 

1.00 

Woods and forests 

.80 

1.00 


^ Adjustments needed for “weeds'* and “grazing.** 

^ For fallow land only, poor condition means “After row 
crop,** and good condition means “After sod,** 


occupied by plant stems or root crowns. Stems and root 
crowns are assumed to reflect the percentage of porosity 
in the top layer tliat is surface-connected by mature 
plant roots to form conduits for air or water. If surface 
continuity is achieved by incorporation or by me- 
chanical means^ the “a** value should be adjusted to 
represent it. 

Estimates of f^ were obtained from the SCS Hand- 
book (^P) wherein major soils of the United States are 
grouped according to their rate of water intake after 
prolonged wetting. Soils are grouped in four hydrologic 
classes: A, B, C, and D. Musgravc(iO) gives the as- 
sociated rates of f^. in inches per hour: A - 0,45 to 0.30; 
B = 0.30 to 0.15; C==0.15 to 0.05; and D = 0.05 to 
0.00. The texture and density of the impeding layer 
give a clue to the selection of f^, within a group, If the 
impeding strata approach clay, f^ is near the lower 
limit of its group; for sand f^ would be near the top. 
Loams would indicate the midpoint fc» 

Depressions on the land surface act as reservoirs 
holding a volume, of water until it is dissipated by 
infiltration, Althougli generally quite small in depth 
(about 0.05-inch for agricultural crops), Vd is drawn 
upon by infiltration and subsequently refilled by rainfall 
for each pause in rainfall intensity. Also, certain agricul- 
tural practices such as level terraces, contoured furrows, 
and listed crops may be designed to give Vd depth of 
several inches, Musgrave and Holtan (31) give concepts 
and some estimates of depression storage in agriculture, 
but generally such, estimates must be based upon 
personal observation^. 
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Certain soils, sucii as the niontmorillonitic days, form 
deep cracks on drying. Crackiiig is estimated for a given 
horizon from the ratio of bulk density at field capacity 
to hulk density when air dry by the equation: 


% Cracks =100 


/ BPW 

\hmj BDD 


[4] 


where 

BDW is bulk density at field capacity in percent 
volume; 

and BDD is bulk density at air dry in percent volume. 

Tills equation was developed by considering a cube of 
soil and by deleting vertical shrinkage. Equation [ 4 ] is 
applied to layer thicknesses to compute maximum 
cracking volumes. The volume of cracks at any ©ven 
time is computed as a linear function of soil moisture 
present, limited to the plant available range, AWC. 
Cracks are therefore at maximum within the root depth 
at wilting point and disappear at field capacity. 

Coefficients for Routing Flows 

Rainfall in excess of infiltration is routed across each 
soil zone of figure 3 and cascaded, subject to further 
infiltration, across designated subsequent soil zones 
enroute to the channel. Overland flow is computed by 
an adaptation of the continuity equation: 


Hoi tan, “ova” can be estimated from vegetative density 
‘'a*’, percent slope “SI”, and feet length “Le” of How 
path by the following equation: 

ova = (150- 125a)'(12/ VIT5) * V^Le [7] 

For laminar flow, “ova” would vary linearly as “SI” in 
equation [7] and the exponent, 1.67, must be changed 
in equation [ 6 ] to 3.0. 

Channel flows and subsurface return flows are routed 
by simultaneous solutions of the continuity equation 
and a storage function. Storage coefficients are obtained 
by integration of the flow recession curve for a given 
watershed. Techniques are under study in the USDA 
Hydrograph Laboratory for the approximation of flow 
recession curves by linear segments to obtain the 
advantages of superposition in routing water flows from 
different areas of the watershed. Flow from each unit is 
routed separately through watershed storage, and all are 
summed to obtain watershed outflow. 

Postramfall recession curves are plotted on semi* 
loganthmic paper in figure 8 for Little Mill Creek 
Experimental Watershed near Coshocton, Ohio. Evapo- 
transpiration must be minimized if we are to obtain a 
true storage-flow relationship. The winter (March) re* 
cession shows the least influence by ET and is used for 
analyses. The equation of a recession curve segment is: 

q, = q^,e-‘/ni [8] 

where: 


[5] 

and 

Qo = {ova)D^“^^ {6] 

where 

Peis the volume of rainfall per At in excess of 
infiltration and depression storage; 

Qq is volume of outflow per At; 

Is overland flow in inches per hour; 
ova is a coefficient dependent upon roughness, length, 
and degree of slope; 

D is tlie average depth of flow in inches; 
and t IS time in hours. 

The coefficient “ova” in equation [ 6 ] is determined 
from runoff recessions on rectangular plots as desenbed 
by Musgrave and Holtan. It is the rate of flow produced 
by an average depth of 1.0 inch of water on the plot 
surface. In the flow data presented by Musgrave and 


Qq is rate of flow at start of a period in inches per 
hour, 

is rate of flow one time-increment later in inches 
per hour, 

m, the absolute value of t/Aln(q), is a constant for 
each straiglit-line segment of the recession curve 
on semilogarithmic papers; 
e is logarithmic base; 

and t is the time-increment in hours. 

By integration of equation [ 8 ] the storage increment » 
AS, within a linear segment of the recession curve is: 

AS = m Aq [ 9 ] 

Values of “m” denved for each linear segment of the 
semilogaritlimic plotting are assumed to represent suc- 
cessive flow regimes (34J, starting with me for channel 
flow and proceeding througli a series-“mi,m 2 , m 3 , and 
m 4 ”-for successively deeper or more devious regimes of 
subsurface flow. Some watersheds, such as the Hastings 
W- 1 1 , have no return flow; hence, only m^ is defined. 
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Figure 8. -Seasonal recession flows, W-97, Coshocton, Ohio. 


At Coshocton, flow regimes must be separated. The 
usual procedure has been to extend segments backmrd 
in time for subtraction from earlier regimes (5). In our 
experience, and tliat of Kill an dais wamy and Secthor- 
man (27), the adjusted points never plotted as a stralglit 
line, and the resulting recession was too rapid. The 
authors prefer to use the same procedure generally 
used for separation overlapping hydrographs of surface 
flows, that is, extending the recession of the first storm 
forward in time for subtraction from subsequent flows. 

Figure 9 illustrates results of tecliniques suggested for 
separating flow regimes for input to this model. Each 
linear segment Js extended toward zero and subtracted 
from subsequent flows to get the rising hydrograph of 
the succeeding regime. Now there is no question about 
maintaining logarithmic linearity vrithin flow regimes. 
Also, the rising hydrograph Fits better with the concept 
that return flows are less when flows in stream channels 
are at higlier stages. The total volume of the separated 
hydrograph can bd computed for successive subsurface 
regimes by the equation: 

S = qb(mL-ni^.i) [10] 


where; 

S Is total storage volume of the difference hydro 
graph in inches; 

is rate at intersect with the previous segment of 
recession; 

m is absolute value of t/ln (qo/qt)i 

L is flow regime number; 

and subscript L -1 indicates of the previous 

segment. 

Input to the model would then be qi> nij, 
rflai q 3 > m 3 , q 4 , and m 4 "’ for those watersheds having 
four discernablc regimes of subsurface flow, The value of 
‘W* is readily determined as the number of hours 
required for the recession segment to cross one log cycle 
divided by 2.3, the natural logarithm of 10 . 0 . Maximum 
free-water storage is assumed equal to the value 
from the soil survey for regimes 1 and 2 , but it is the 
volume under the hydrograph computed by equation 
[ 10 ] for the remaining flow regimes as in figures 8 and 
9 . 

M^,,as determined from the watershed flow recessions, 
includes detentions of overland flow. Since overland flow 
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Figure 9. -Separation of flow regimes, W«'9?, Coshocton, Ohio. 


is routed in equation [6], the depth of overland flow 
must be subtracted from the AS of equation [9] . At the 
peak, outflow is essentially equal to inflow arid channel 
storage would be. 

S = * q -(q/ova) [11] 

We can get an approximation of the corrected 
(cm^ in the program) by solving for m in equation [9] at 
one unit of flow rate. 

cm^ rn^. - ova*"'^ [12] 

Equation [12] is an approximation only and typifies 
an area of needed research. 

Concern for the transport of agricultural chemicals 
has placed emphasis on the path of both surface a?id 
subsurface flows from relatively small upland areas well 
above the elevation of emerging subsurface return flow. 
As a result* values of mv to are obtained from 
downstream gaging sites of an encompassing watershed 
together with an indication of the number of regimes 
(including 111 ^) which occur above the weir. Output then 
includes ‘‘onsite’* return flow in “runoff’ and lists 
“offsite” flows for comparison with regional informa- 
tion. 


As the size of watershed diminishes, the importance 
of detail increases. Therefore, mi and are subscripted 
by zone and reevaluated from soils information. As- 
sumptions were: 

(1) Outflow is equal to permeability - width/length 
• gradient; 

(2) permeability varies as the ratio of pore diameters 
squared or as the square of the ratio of large pore 
volumes; 

(3) and, since is permeability of layer 2> the 
permeability of layer 1 is fcCGi/Ga)^. 

Then, since storage, (FREE WATER) is equal to 
DEPTH/G and since m«= FREE/G, we can derive 
and mz(2) for each zone by equation [13] . 

mz(L) = depth (L) • G 2 ^ • Le^ • (Sp + 

(f^ • G(L) • zone area • SI) [13] 

If tlie downstream flow record indicates that the 
hydrograph recession passes abruptly from the character- 
istic surface flow hydrograph to the low, sustained 
baseflow without transitional slope breaks, mi and 
can be entered as zeroes. Tlie program will not compute 
flows in a regime whose m value is zero. 
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Equation [14] is a combination of the continuity 
equation and equation [9] used m computing outflow 
from each regime: 


2AI 

” 2m + At 


2m “ At 


{14] 


Equation [9], being linear, forces the routed hydrograph 
from equation [14] to peak at its crossing with the 
inflow hydrograph, whereas channel flow responds non- 
linear ly to changes in inflow rate. Therefore, non- 
linearity was induced on the rising side of the 
streamflow routing by the techniques of Holtan and 
Overton (21) for successive routings through one-half of 
indicated storage using equation [15] : 


, 2AI f m - At 

” m + At ^ m + At 


[151 


In the second application of equation [15] the quantity 
Cqi V q? OAt/2 is substituted for AL The symbols are: 


I is inflow volume in inches; 
m is routing coefficient in hours; 
t is increment of time in hours; 
q is rate of outflow from storage in inches per hour; 
is theoretical rate of outflow from first half of 
storage in inches per hour; 

and subscripts i and 2 indicate beginning and end of 
At. 


Ground water recharge from the ultimate return-flow 
regime is estimated on a regional basis Average annual 
rainfall, average annual ET, and average annual stream- 
flow yields in the vicinity can be used to derive an 
average annual ground water recharge. TMs can be 
converted to inches per hour as an estimate of gr. At 
Hastings, this value was estimated as equal to but the 
estimate for the Coshocton, Ohio, watershed was 0.0006 
inches per hour, and for the Texas and Florida water- 
sheds, gr was estimated as zero. 

It must be borne in mind that for small upland areas 
the deeper regimes may pass beneath the point of 
surface-flow measurement. If so, gr must be estimated 
equal to of the first unmeasured regime and the q 
and m values of the nonoccurring regimes must be 
omitted. 

Increments of downward seepage, ‘‘subput'', to the 
next regime, are computed as a function of a free-water 
present: 


SubputL+i = t - C ^ (G “ SA)/G [17] 

The potential rate of outflow qp from a regime is 
derived by equation [14], but actual outflow from a 
regime is held equal to zero, as illustrated in figure 9, 
whenever summation outflow exceeds the value of 
intercept for that regime. When summation outflow 
is less than a favorable gradient is assumed and 
outflow q|^ is computed by the following equation: 


Hydrogeology 

In the USDA Hydrogiaph Laboratory’s model, in- 
filtrated water is proportioned to ET, or free-water 
evaporation, and to downward seepage, or to lateral 
return flow, in each flow regime. Downward seepage and 
lateral flow are supplied by free water; hence, estimates 
of maximum seepage rate “C” and free water capacity 
as well as the storage coefficient “m”, are needed 
In each flow regime. 

Since subsurface flow regimes are considered se- 
quential, seepage, C, from a given regime is inflow to the 
next regime and must bne adequate to supply the sum of 
the maximum flow rale in the next 

regime and the rate of deep seepage, gr: 

where: 

L is the regime number; 

and gr is the maximum rate of ground water recharge 
in inches per hour. 


Al flp (fib ^q)/fib {1®] 

where 

qjj is intercept mlh previous segment; 
is sum of all existing flow rates; 

and qp is potential regime outflow from equation 
[14]. 

An exception to equation [18] is made for frozen 
ground. If the 3-week moving average air temperature is 
below freezing, the c of equation [17] is reduced for 
layer I by 0.02 inch per hour for each average degree 
Fahrenheit below 32. Thus it is assumed that liighly 
permeable soils become impervious at lower tempera- 
tures. Since only fc is affected, soils take in water up to 
saturation during low temperatures but no outflow 
occurs from layer 1 until temperatures rise above 
freezing. 

The limits of free-water storage in the first two 
regimes are computed from soil survey data as the 
products of horizon thickness and percent freely drained 
porosity. The limit of free-water storage, G, associated 
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with subsequent flow regimes is computed from tlio 
storage^flovv relationship: 

Gl =qb * 119] 

There is no assurance that the flow recession analyzed 
represents a saturated A or B horizon. Hence, the 
G-values determined from soil survey data for these twc 
horizons lake precedence in the model. Subsequently, 
the q-values and (he G-values of regimes 3 and 4 are 
apportioned to each zone of ilie watershed in proportion 
to the frec-wa ter limit “G’^of the topsoil: 

Ga (zone) = mj • ' Gj Czone)/Gj [ 20 ] 

where is a vveiglitcd watershed average from soil 
survey. 

By (his means the storage indicated by the lower flow 
recession is dispersed over the waterslicd in proportion 
to the depth of the topsoil horizon m each zone or 
hydrologic response unit. Routing coefficients, nij and 
1114 , remain unchanged for each regime, but the storage 
potentials, and the maximum outflows, are 
now unique for each zone. 

qmav = G/ni [21] 

Evapotranspirallon is drawn from both G and AWC 
wlien roots are present. Geologic research is needed to 
develop concepts of flow regimes associated with spe- 
cific depth in the profile. In the meantime, we are 
placing greatest leliancc on soil surveys for storages of 
regimes 1 and 2 and depend on the recession curve 
analyses for storage values of subsequent regimes. Roots 
are arbitrarily designated as drawing water from the first 
two regimes, often creating a deficit In AWC to be 
satisfied before accruals to subsequent regimes can 
occur. Therefore, A\VC| and AWC 2 are required input. 
AWC I IS obtained from the moisture tension data of the 

I 

first layer. AWC^ is estimated from the moisture tension 
data below the first layer to the depth of suitable habitat 
for roots as indicated in the soil profile description. 
Moisture deficits are computed only in the depth 
indicated for root habitat. Today, there is no basis for 
distributing any part of AWC to subsequent regimes 
except in a few places like Florida where records of wells 
over impervious rock at very shallow depths may permit 
some speculations. 

Evaporation from free water in the soil is treated as a 
function separate from evapotranspiration. It is ex- 
tracted from the sum of free water “G” in regimes 1 and 
2 as the product of percent total free water available and 
pan evaporation. 


A percentage of overland flow tioin each zone is 
designated to cascade across the subsequent soil zone, 
with the remainder, if any, allocated to the alluviums or 
directly to channel inflow Allocation is a matter of 
judginent based on the nearness of a zone to aJliiviums 
or upon the presence or absence of diversions, graded 
terraces, listed furrows, or other interceptors leading 
directly to a stream channel. Infiltration in excess of 
AWC in the root zone is routed through storage of 
subsequent regimes and accumulated as run-on to the 
alluvium or as direct inflow to the stream channel if so 
designated by input. The total of all inflows to Uic 
channel is routed, using the channel storage coefficient 
in equations [13] and [15] to obtain the outflow 
hydrograph for the watershed. 

The Erosion-Deposition Submodel 

'"A Watershed Sediment Yield Model with 
Hydrological and Chemical Interfaces" 

This submodel predicts components necessary for the 
erosion-chemistry interface. It simulates the erosion- 
deposition continuum for watershed slopes on an area 
basis. This phase is subdivided to estimate the interrill 
and rill erosion components of the total sediment load. 
Particle size distributions of the eroded material are also 
estimated for erosion-chemistry relations. Each calcula* 
tion utilizes hydrologic input together with basic soil 
parameters for each compartment of the watershed 
transformation. 

Basic Equation 

The conceptual basis for the erosion and sedimenta- 
tion submodel is the Universal Soil Loss Equation 
(USLE): 

A=RKCPSL [22] 

where 

A == estimated soil loss (T/ac) 

R = rainfall factor 
K = soil erodibility factor 
C = cropping management factor 
P = erosion control practice factor 
S = slope steepness factor 
L - slope length factor 

The S and L factors can each be expressed mathe- 
matically based on data analyses by Wischmeier and 
Smith ( 55 ), 
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^ 0.43 + 0.30s +0,043s^ 

[23] 

^ 6.613 

XO.5 

L = (X/72.6)0 5=A_ 

[24] 


where 

s percent slope 

X” slope length 

The USLE is widely used, and each of the factors has 
been determined for many soils and situations throu^v 
out the United States. It was originally derived to 
predict average annual soil losses. However, all factors 
can also be evaluated on a storm basis. Throughout this 
report, it is being used on a storm basis. In the USLE, 
the energy available for detachment is that derived only 
from raindrop impact. 

Williams (5.?) has modified the USLE to piedict 
sediment yield, G, at the watershed outlet. The rainfall 
energy term was replaced by an empirically derived term 
based only on flow characteristics for a particular storm 
event. 


neously at times during a rainstorm. The relative 
influence of each is dependent on both time and space in 
a watershed. 

At the beginning of a storm, before runoff begins, 
raindrop splash is detaching and transporting soil. This is 
the major energy source since runoff has not begun. 
When runoff begins, it also has potential that increases 
with the rising hydrograph. At the end of the storm, 
runoff is the detaching and transporting agent. 

With regard to position on a watershed, upslope 
erosion, where detachment is limiting, depends primarily 
on rainfall energy and slope gradient {54), Splash action 
and sheet flow transport the soil particles to flow 
concentrations or nils. Downsiope, as runoff increases, 
its detachment and transpoit potential increases relative 
to rainfall. Finally, transport capacity becomes limiting 
and sediment yield becomes closely related to the 
hydraulic properties of the runoff and less related to 
rainfall characteristics {Wa\ 

Foster (ida) modified equation [1] to reflect both 
rainfall and runoff contributions. 

A = [aR + bcQqpl/3] kCPSL [26] 


G = 


a(Qqp)^ 

~Tr 


KCPSL 


[25] 


wliere 

G = sediment yield (T/ac) 
Q - runoff volume (ac-ft) 
qp = peak flow rate (efs) 

A^ = drainage area (ac) 
a, P = coefficients 

K, C, P, S,L- USLE factors 


The coefficients, a and ft determined from 778 events 
on 18 watersheds in Texas and Nebraska, were 95 and 
0.56 respectively {53), The coefficients, calculated from 
data by Piest, et al. {37) from 103 events on two Iowa 
watersheds were found to be 23.6 and 0.65, respectively. 
The exponents are nearly equal, indicating that a change 
in the product, Qqp, has about the same effect. 
However, location differences are revealed more isignif- 
icantly by differences in a. 

The USLE and the Williams modification represent 
the two extremes with regard to energy required to 
initiate detachment, in that one utilizes rainfall energy 
entirely and the other, runoff entirely, Both rainfall and 
runoff have detachment potential and act simulta- 


where 


a,b = weighting parameters (a+b= 1) 
c = equality coefficient 
R = rainfall factor 
Q == runoff vol. (in.) 

Qp = runoff rate (In, /hr.) 

K, C, P, S, L = USLE factors 


The weighting parameters reflect the relative amounts of 
erosion caused by rainfall and runoff under unit condi- 
tions (L, S, C, P= 1). Limited experimental evidence 
{56y 2S) indicates that the detachment is about evenly 
divided under these conditions (i.e,, a=b=,5). The 
equality coefficient, c, can be estimated from rainulator 
plot data under unit conditions. 


c = (2AjK-R)/Qqpl/3 [27] 


where 


- soil loss under unit condition. 

Analysis of 55 unit rainulator runs on 20 soils in 
Indiana*^ and Minnesota shows that c equals 29.3 or 


"^Foster, G. R, 1973, Personal communication. 
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approximately 30. Substituting these values of a> b* c 
into equation [2b] yields the basic equation for the 
erosion-deposition submodel. 

A = EKCPSL [28] 

where 

E = 0.5R + ISQqp*/^ 

The coefficients, a, b, and c, may be dynamic in that 
they may change from time to time and place to place 
on a watershed during a storm. They may also be 
different during the detachment and transport phases of 
erosion. So might the other factors in the USLE, 
However, formulations of these relations- cannot be 
made at tliis time due to lack of data. Therefore, 
throughout the succeeding discussion, it is assumed that 
these coefficients and the USLE parameter remain 
constant at their indicated values througlt time and space 
and also for both phases of the erosion process. 


Detachment Component Development 

For a field of uniform slope, s, as shown in figure 10, 
equations [24] and [28] are combined to give: 


Foster and Wischmeier {10) have shown tlv.it the 
detachment rate at any point x is 


D 


* 3x 


or from equation [30] 

EKCPS (1.5 x°-^) 

185.58 


[ 31 ] 


[32] 


Integration of equation [32] between two length .limits 
yields tlie detached soil capacity for the lengtli segment. 
So, for segment j of the complex slope shown in figtire 
11 


Adj = dx [33] 

where = detached soil capacity for segment j or 

185.58 ' ^ 


where 

E..., * 'SQ.,., (W,)’" 


EKCPS 

8.521 


When this is extended to a compound slope having n 
segments, the cumulative detached soil capacity is 


This can be extended to any point x on the slope so that 

. EKCPS x®'^ [29] 

^ 8.521 

When the soil loss is expressed as Ibs/foot of width, 
equation [29] becomes 


At = Adi + • • • + Adj_i + Adj + Atij+ 1 + • • • + Adn 

This would equal the sediment yield per unit width if no 
deposition occurred along the slope length. Alterna- 
tively, tlie sediment yield equals the cumulative de- 
tached soil capacity provided the transport capacity is 
not limiting. 


EKCPS 

185.58 


[30] 
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Transport Component Development 


For a uniform slope» let the USLE represent the 
transport capacity at any point x, so that 


A 


cx 


CKI'S 

185.58 


[35] 


where 


= transport capacity at x 


This assumes tliat C, P, and S remain constant and 
identical to those for detachment in the transport phase. 
The value of K should reflect tlie transportability of the 
material coming from ups! ope and tlie effect of soil 
properties on rilling and transport capacity on down- 
slope segments. Tlrerefore, an effective value, K, is used 
which is the average K of the upslope segments. For a 
nonuniform slope, the transport capacity for any seg- 
ment is assumed to be that sediment load calculated 
using the total slope length for L and the slope steepness 
of the segment for calculation, S, Refeiring to figure 11, 
tlie transport capacity at Xj+j is 


185.58 


[36] 


assuming a slope steepness of Sj througliout the lengtli 

Therefore, two cases can exist for slope segment j 
(i,e,, at xj+i). 


Case 1 

If A^<A^,, erosion occurs on slope segment j 
because the transport capacity exceeds tlie total de- 
tached soil, Therefore, the sediment yield 
G = A^i + • • • A^^j and the deposition M == 0. 


example, a chemical application tilled-in would have a 
different contribution from another applied as a spray. 

The conceptual basis for the nlMnterrill contributions 
are separate laboratory studies conducted by Young (5d) 
and Meyer, et al. {28) in which they observed the trends 
illustrated in figure 12. The top part of the figure shows 
that the interrill erosion rate builds up rapidly and then 
abruptly levels off, indicating that the rate remains 
constant downslope. The bottom part indicates that at 
some distance downslope rill erosion volume equals 
interrill erosion volume and, again, that interrill erosion 
increases linearly. Furthermore, when interrill and rill 
erosion volumes are equal, interrill erosion rate has 
already become constant. 

Foster*^ derived a relationsliip which approximates 
the ratio of lill erosion to interrill erosion fora segment 
where detachment rate is at capacity. 

(x|-x?_i) (0.043s?) 

“ 72.6(Xj~Xj_i ) (0.3sj +0.43) 

[37] 



where s = percent slope and 
Kj./K| = ratio of rill to interrill soil erodibility 
(raiibility). 

Rough guidelines for the ratio, are: 

Kj/Kj = 2; “severe rilling” 

1; “moderate rilling” 

Kj/Kj = 0.5; ‘little evidence of rilling” 

In general, 0 increases as slope length, steepness^ rllli- 
bility and runoff increase. 

The quantity of sediment detached by rill erosion on 
a segment is 

[38] 


Case 2 


and the quantity detached by interrill erosion on a 
segment is 


If deposition occurs on slope segment j 

because total detachment exceeds the transport ca- 
pacity. Therefore, the sediment yield G = 
deposition M A^ - ^cxi+j * 


lAj-Ajj/Oij^l) [39] 

Tlie equations for cumulative rill erosion and Interrill 
erosion in a compartment are: 


Rilt-Interrill Contributions 


m ni 

RA= S RAj andIA= 2 lAj [40] 

j=i j=i 


The division of soil loss between rill and interrill areas 
is important for the chemical-sediment interface because 
of different methods of chemical application. For 


19 


where m is the number of slope segments. 



[42] 





RD = 


324 (RA) 
wX 


where RD = rill depth (in.) 

Eroded Material Particle Size Distribution 

It is important that not only tlie total amount of 
sediment be predicted, but also tJic composition of 
eroded mateiial should be known for the sediment- 
chemistry interface. This problem was approached utiliz- 
ing texture information from 56 Midwest soils. The 
original soil textures ranged from sand to silt loam to 
clay. The data consisted of measured soil loss and 
textural classification of both the soil matrix and the 
eroded material. This was done for each of two 
rainulator runs on four replications for each soil. 

For purposes of analysis, the following specific 
surfaces were assumed (i7); 

clay “ 200m^/g 


DISTANCE OOWNSLOPE (ft) 


silt - 40m^/g 


Figure 1 2.~Eroded soil apportioned between nil and interdil 
sources for nil-susceptible Plot R-3 (From Meyer et. al. (28),) 


Tlie interface to the chemical submodel requires that 
rill surface area and depth be specified. Limited experi- 
mentaJ evidence obtained by Foster^ on rainulator plots 
indicates that the percentage of the surface intercepted 
by nils at vanous cross-slope transects remains relatively 
constant from the top to' the bottom of the slope, 
although the number of nils becomes less in the down- 
slope direction. Based on these data and other data 
from construction site erosion, he postulates that the 
average liO width can be represented as follows; 



[41] 


sand - 0.5ni*/g 

Tlierefore, the specific surface for a particular soil, i is 
SSj = 200(%C1) +40(%Si) + 0.5 (%Sa) [43 J 


The clay ratio can also be calculated by 
%C1 

" %Si + %Sa 


[44] 


Figures 13 and 14 show the graphs of the specific 
-urface and clay ratios of the soil matrix versus those for 
the eroded materials for the 56 Midwest soils. The 
regression equations are as follows: 


SSe == 14.6 + 0.84 SSm 


[45] 


where 


and 


w * average rill width (in./ft.) 

This indicates that the liii width depends on the 
rililbility of the soil, Kj/Kr, and the slope steepness, s. As 
rillibility and slope steepness increase, the average rill 
width increases. 

Average rill depth can be calculated by using the 
results of’equaUons {40] and [41] and utilizing the fact 
that 1 yard^ of soil weighs about 1 ton assuming the 
bulk density of 1.2, 


CRe = 0.021 + 1.08 CR;n [46] 

where the subscripts m and e refer to soil matrix and 
eroded material, respectively. The coefficients' of de- 
termination for equations [45] and [46] are 0,81 and 
0.82, respectively. Equations [43] -[45], together witJi 
the fact that the sum of the components equals 1, 
constitute the estimation procedure. This is a first 
approximation recogni 2 dng that many factors are not 
included. It is illustrated by the following example. 
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Example: Estimate the mechanical composition of the 
eroded soil if the composition of the soil 
matrix is 

“ 20% “ 40% — 40% 


CRni 


0,20 

0,40 + 0,40 


0.25 


From equations [45] and [46] the estimated parameters 
for the eroded materials are: 


From equation [43] 

SSm = 200(0.20) + 40(0.40) + 0.5 (0.40) = 56 
and equation [44] 



SOIL MATRIX 

Figure 13. -Plot of specific surface ; eroded vs. soil matrix. 


SSc = 14.6 + 0.84(56) = 62 

CRe = 0.021 + 1,08(0.25) = 0.29 

Therefoie, solving the following three equations for the 
eroded soil textural components 

200 Clc + 40 Sic + 0.5 Sae == 62 

= 0.29 

Sic Sac 

Cle + Sie + Sac " 1 

sliow that the estimated eroded soil compositions are 

Clc = 22% 

Sic = 43% 

Sae = 35% 

and the enrichment coefficient is 
Cle 22 

147] 



CR^ — SOIL MATRIX 

Figure 14. -Plot of clay ratios; eroded vs‘. soil matrix. 


Watershed Calculations 

The orosion-depositlon submodel is applied to the 
waterslied utilizing the transformation illustrated in 
figure 2. Each compartment may have a unique set of 
parameters required for the solution of the 
erosion-deposition equations. Calculations proceed from 
tlie top, downslope to tlie channel, continuing from one 
zone to the next. It is assumed that all soil delivered to 
the channel reaches the watershed outlet. 

The hydrologic submodel predicts tlie flow volumes 
and discharge rates from each zone but not from cacli 
compartment, Apportionment of flow volumes and peak, 
discharge rates among the compartments within a zone 
were made using linear interpolation. The same tech- 
nique was used for adjacent compartments between 
zones. Because the zone is the smallest calculation unit 
permitted by the hydrologic submodel, all hydrologic 
input variables are Identical for compartments traversing 
the slope. Input variable differences along the slope can 
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be handled by adding more zones to the hydrologic 
niiodeh 

The Chemical Submodel for ACTMO 

''A Model of Chemical-Soil-Water Interactions” 

The objectives of tlie chemical submodel arc to 
calculate concentrations of a chemical in runoff, the 
total amount of the chemical leaving the watershed, and 
tile location of the remaining chemical, so that manage- 
ment practices can be selected to minimize undesirable 
movement, The submodel traces the movement of a 
single application of a chemical tlirough and over the 
watershed. Cultivation is the only management practice 
explicitly accounted for since it redistributes the chemi- 
cal. All other practices such as contouring and terracing 
are implicitly handled tlirougli the hydrology and 
erosion submodels. Nitrate behavior was judged to be 
unique enough to warrant a separate submodel that uses 
many of the features in the general model, 

Simple chromatographic theory is used to describe 
tiie leaching of the chemical and to provide a method of 
calculating surface concentrations for runoff and ero- 
sion. Cliemicals are degraded or lost, redistributed by 
cultivation, and nitrate is mineralized and taken up by 
the crop during the period of time between storms. Dur- 
ing the storm period the chemicals are transported by 
water and sediment. Adsorption characteristics of chem- 
icals range from highly adsorbed to no adsorption. Highly 
adsorbed chemicals, such as DDT and phospliate, have 
very low solution concentrations and are moved princh 
pally by eroding the soil on which tlxey are adsorbed. 
Other chemicals, such as nitrate, are not adsorbed and 
move with the water. Chemicals with intermediate 
degrees of adsorption are moved by both vehicles. 

The multitude of chemical reactions in the watershed 
are classified into three groups: adsorption or the 
distribution of the chemical between the solid and the 
solution phases, breakdown or loss of chemical, and the 
dispersion of the chemical as it moves through the soil, 
Each of these three groups of reactions is described by a 
simple relation which in many cases is only a crude 
approximation. More complete descriptions require in- 
formation not generally available. 

The compartment is the basic unit for computations. 
The hydrology submodel divides the watershed into 
ones while the erosion submodel creates routing tubes 
carry the cascading water and sediment across the 
to the channel. The compartment is tlie soil 
V of each zone within a routing tube shown in 
of the watershed transformation section. For 
chemical transport, the hydrology submodel 


is limited to a single ciop per zone and assumes the soils 
are similar within a zone. Althougli not essential, the 
layers of chemical interaction with the soil are assumed 
to be tlxe same as the hydrologic layers. Tlie English 
units, such as pound, acre, and inch, are used. A 
bibliography of chemical interaction and chemical trans- 
port (Appendix) has been compiled during the develop- 
ment of this submodel. These papers should be of 
assistance to anyone wishing to estimate the parameters 
for a given chemical and soil or tliose wishing to modify 
any concept used in this model. 

The Chemical Option 

Adsorption.— The adsorption relation is assumed to be a 
linear isodierm, 

S = AC • C [48] 

where S is the pounds of chemical adsorbed per pound 
of solid, C is the pounds of chemical per pound of 
solution, and AC is the adsorption coefficient. Batch 
equilibrations and column leaching studies^( Appendix) 
provide information on the magnitude of the adsorption 
coefficient for many soils. The model can be easily 
changed to include nonlinear relations, but this requires 
another parameter. 

Rates of chemical adsorption and desorption {26) 
were considered to introduce more complexity than is 
merited for this simple model. Nonequilibrium condi- 
tions cause tlie chemical distributions to be asymmet- 
rical, and are probably more prevalent in tlie real world 
than not. However, it would appear that this effect is 
usually not c'lny larger tlian the normal variability 
encountered in the field {40, 41), 

Breakdown.— There are numerous papers (Appendix) on 
the breakdown, degradation, loss by volatility, and 
persistence of chemicals in both laboratory and field 
conditions. A simple first order rate equation was chosen 
to represent these many processes, The chemical remain- 
ing, A #/a, is related to the initial amount, AA #/a, by 
an exponential function (EXP) of time, TI days, 

A = AA • EXP(-BC • TI) [49] 

The rate coefficient, BC/day, is dependent upon tlie 
environment around the chemical (i). There isn’t 
enougli quantitative information available at this time to 
make tlxe necessary corrections, but tlxe model has been 
designed to accommodate the calculations when the 
relations are developed. The average temperature and 
soil moisture are available from the hydrology submodel. 


22 



Tlie type of application, banding or broadcasting, also 
changes the degradation rate (d) as can repeated applica- 
tions of the same chemical (id). Such factors have been 
neglected in this first approxirnation.Xarge amounts of 
pesticides can be lost by volatilization during application 
and, iheiefore, the amount at the start of the simulation 
is not tlie amount leaving the applicator, but rather the 
amount that actually reaches the soil. Volatilization loss 
after the chemical reaches the soil is usually significant 
only if othei losses are small. Volatilization is a complex 
process involving a number of factors, and this first 
version assumes that these losses will be reflected in the 
apparent breakdown coefficient. A detailed discussion of 
pesticide volatilization is given by Spencer et al, {43\ 
Dispersion.-“When a chemical is carried through the soil 
by flowing water, there is a broadening of the chemical 
band and a decrease in the peak concentration. Two 
processes are in operation: diffusion and dispersion. 
Diffusion is the movement of the chemical from areas of 
higli concentration to areas of low concentration. 
Hydrodynamic dispersion is the lesult of unequal water 
flow through the different sized pores and cracks in tlie 
soil and should reflect the aggregate size {36) and the 
ratio of porosity to permeability {4). The dispersion 
coefficient has the same units as the diffusion coefficient 
(area/time) and is proportional to the waterflow 
velocity. For tlie flow rates during infiltration and 
drainage, the dispersion coefficient is considerably larger 
tlian tlie diffusion coefficient. A more useful factor for 
tlie equations to be used is the ratio of the dispersion 
coefficients and the average flow velocity. For the sake 
of clarity this is called the dispersion distribution 
coefficient, DC, with units of length. 

Equations involving dispersion are matliematically 
similar to diffusion equations where many solutions are 
available (24). Some experimentally measured values of 
DC from column studies range from .03 inch with 
sieved soil (7) to about 5 inches with undisturbed soil 
columns (7^). Little field information is available for 
tins parameter. 

Movement.— There are basically two ways to calculate 
file movement of a chemical through the soil. Tlie first is 
the finite element method in- which the soil is divided 
into thin layers of soil (9, 1 2). The second is to assume a 
simple continuous chromatographic process. In the first 
method the chemical is moved from one layer to another 
and complete mixing with uniform distribution within a 
layer is assumed, Two problems arise with this approach. 
Wlien tlie chemical is applied in a narrow band, a large 
number of layers are needed to accurately trace tlie 
chemical movement. The second problem is to define 
the thickness of the soil surface-runoff interface. The 


ratio of the watershed surface area to the specific surface 
area (area/area per unit weight of soil) is one estimate of 
tlie amount of soil in this interface (25), but this appears 
to be a considerable underestimate. 

An effort was made to circumvent these problems by 
assuming an idealized bell-shaped distribution or plug 
distribution for use in a simple chromatographic process. 
This assumption of a known mathematical equation 
provides a means to calculate the concent lation at any 
deptli including the surface of the soil. It appears that 
such a procedure is adequate when there is no adsorp- 
tion or the adsorption isotherm is linear. This holds for 
nitrate, phosphate, and many pesticides, but not for 
cation exchange involved in salinity movement (39). 

Water moves the chemical tlirougli the soil along 
several patlis. The infiltrating water tends to move the 
chemical vertically and sometimes laterally througli the 
soil. Runoff water carries some cliemical with it and also 
causes erosion which can carry additional amounts of 
chemicals. Table 6 gives an outlme of the calculations 
made in linking chemical characteristics with the move- 
ment of tlic carrier. Two types of chemical applications 
are provided for on the surface or mixed into the soil. 
Tire amount applied is assumed to be corrected for 
application losses. The surface application includes the 
amount on the plant surface, if any, that is readily 
washed off. A major gap exists in the understanding of 
factors controlling the rate of washoff from the plant 
surfaces. Mixing the chemical into the surface soil is 
assumed to result in an ideal uniform plug-type distribu- 
tion even tliough this may not always be true. 

Table 6.-Outline of steps for tracing chemical 
movement 

1. Amount of loss before the storm. 

2. Movement by infiltration, 

3. Concentration at the surface. 

4. Runoff, 

5. Intcrrill erosion. 

6. Rill erosion. 

7. Accumulation and transport of water, sediment, and 
chemical along routing tubes. 

8. Deposition, 

9. Lateral outflow. 

10. Mixing of water, sediment, and chemical from other 
tubes, 

1 1. Location of the remaining chemical 

Losses before storms, step 1.— The chemical is lost 
according to equation [49] for tlie period of time from 
application until the first storm and between storms 
thereafter. The time-average soil moisture and tempera- 
ture are available for tliis period and could be used to 
modify the amount of breakdown or loss. An extra 
coefficient for the pliotodecomposition rates at the 
surface is used in the array of breakdown coefficients. 
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Infiltration, step 2.—A great deal of information has 
been developed about the infiltration process (/J). Since 
it was decided to make this first attempt a very simple 
inodelj the leaching process was divided into two stages. 
The first stage is prior to runoff when nonadsorbed 
chemicals like nitrate and halides move with the water 
front into dry soil (55), This distance is approximately 
die amount of infiltration in inches, IN, divided by the 
volumetric field capacity, FC. 

D = IN/(FC+BD-AC) [50] 

The term (BD*AC) accounts for the chemicals lagging 
behind the waterfront when the chemical is distributed 
between the water and the soil of bulk density BD ( 13 ), 
Some anions like nitrate are excluded from some of the 
fine pores in the soil, and BD‘AC can be given a negative 
value corresponding to the water volume from which the 
anions are excluded ( 42 ). The term field capacity is used 
as an operational definition of the amount of water 
remaining after saturated flow has ceased. It Is recog- 
nized tliat tliis is physically imprecise and difficult to 
evaluate. 

In this model the second stage of leaching is for the 
infiltration that occurs during runoff, Tltere is consider- 
able data showing that continuous percolation is less 
efficient in leaching salts than intermittent flow. Tliis 
inefficiency is probably a result of relatively more water 
moving through large pores than small ones. An ex- 
amination of the data in Miller et al, ( 29 ) shows that the 
inefficiency is proportional to the water content above 
field capacity. Therefore, tire leaching efficiency used is 
the ratio of field capacity to the average moisture 
content during this period. Tlie average moisture content 
is calculated as halfway between field capacity and 
saturation since the process of reaching saturation and 
the return to field capacity accounts for the bulk of the 
water movement. The depth of chemical movement, D 
inches, is calculated by 



FC 

AM+BD-AC 


[51] 


1/4 in/hr. When the soil was air dry, the bromide, which 
is leached similai to nitrate, moved with the water front 
and was found at a depth of 4 inches in a sharp band. 
Wlien the soil was wet at nearly field capacity, the 
simulated rain soon ponded and suction was used to 
draw the water througli the soil. In this case there is a 
loss of leaching efficiency and a more dispersed band is 
found at a depth of 2.5 in. Lithium is a sliglitly adsorbed 
chemical with an adsorption coefficient of about three 
in tliese expenmental conditions and thus is moved only 
about one-tenth of the bromide distance. Tills distance is 
too small to show any effect of the wet or dry 
treatment. 

The concept being used for leaching efficiency and 
dispersion was tested for field conditions with publidred 
data of Miller et ah ( 29 ), Potassium chloride was applied 
to small plots and leached using three irrigation treat* 
ments: continuous ponding, successive 6-inch incre- 
ments, and successive 2-inch increments. The chloride 
content of the water at different depths was monitored 
with porous cups and soil water by tensiometers. The 
dispersion distribution coefficient was calculated from 
all the peak concentrations relative to the distance they 
had moved. The calculated peak location, height, and 
width are compared with five sets of observed data in 
figure 16. The fit with othei observations is equally good 
which suggests that the relations used should be ade- 
quate for the first approximation model being de- 
veloped. 
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where IR is the inches of infiltration during runoff AM 
is the average moisture content, and the other variables 
have been defined. 

Experimental evidence has been obtained to verify 
some of these relations. Figure 15 illustrates the be- 
havior of surface-applied salts during infiltration on soils 
')f different moisture contents (unpublished data from 
)urant, Oklahoma), LitJoium bromide salt was applied to 
le surface of columns of a Fort Collins silt loam soil 
then sprinkled with one inch of water at the rate of 


4 - Br tdryj 

- - i 1 I I 1 I 

0 12 5 

ME/ 100 a SOIL 

Figure 15. -The movement of lithium bromide through a column 
of Ft. Collins silt loam by one inch of simulated rain when 
the antecedent moisture was air dry or at field capacity. 
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Concentration at different depths and times, step 
3.-The combination of simple chromatographic be- 
havior with dispersion and a Gaussian or bell-shaped 
distribution gives rise to equation [52] from which the 
solution concentration, C ppm. at any depth (X inches) 
can be calculated. 


C(x) = 


A-U 

\/4rr*D(>D 



[52] 


A is the amount of chemical and U is a units factor 
equal to 4.42/(SM+BD* AC) when A is in pounds per 
acre and C is in ppm at a volumetric moisture content of 
SM- When X = 0, the equation gives the concentration 
at the surface and for X = D the peak concentration in 
the soil is calculated. 

Wlien the chemical is mixed into the soil AP inches 
deep and there is infiltration that moves the plug 
downward D inches, the boundaries of the plug are 
dispersed in a sigmoid pattern (55). The concentration at 
X inches from the middle of the plug is given by (24): 



where ERF is the error function and Cl is the initial 
uniform concentration in ppm: 


Cl - A/(.3333 AC + ,226^ SM)AP [54] 

with .3333 million pounds of soil per acre inch and 
,2266 million pounds of water per acre inch. Equation 
[53] is used to calculate the concentration at the 
surface, X - D + 1/2 AP, and to check if the plug 
distribution can be treated as a bell-shaped distribution. 
When the concentration at the middle of the plug X = 0 
has been depleted more than five percent by the 
dispersion effects, then the distribution is assumed to be 
bell shaped, 

After a number of storms the concentrations in the 
solution and on the solid at the surface may be reduced 
below a predetermined critical level, and the runoff and 
erosion calculations can be bypassed. 


Runoff, step 4.— Affer the initial infiltration when the 
water begins to pond on the surface and runoff 
commences, it Is assumed that the rainfall ^equilibrates 


PARTS PER THOUSAND CHLORIDE 
O I 2 3 4 50 I 2 3 



Figure 16. -The calculated and observed solution conconlrntlons 
of chloride with depth under two irrigation treatments 
(Panoche clay loam, California), 


with the solid phase on the soil surface. Then the 
concentration In the water as it moves away from the 
surface is the same whether it moves vertically as 
infiltralion or laterally as runoff. Under these conditions 
a projected depth of Icacliing can be calculated with 
equation [51] using the sum of infiltration and runoff. 
The concentration in the solution at the surface, and 
thus in the runoff, can be calculated with the projected 
depth In equations [52] or [53] . 

The loss of the chemical to runoff should be 
calculated by integrating the surface concentration as 
the depth increases to the projected depth. No simple 
way has been found to perform this integration so an 
average concentration is calculated from equations [52] 
or [53] using an average distance between the depth at 
the start of runoff and the final projected depth, 
Examination of the behavior of these equations indicates 
tliat there will be an underestimation of the mtan 
concentration when the runoff and infiltration remain 
constant. However, the infiltration is usually higher at 
the beginning of the storm, carrying a larger fraction of 
the higli concentration into the soil, which reduces the 
underestimation. The total amount of chemical lost in 
runoff Is calculated as the product of the mean 
concentration during the runoff and the amount of 
runoff. 
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interfill erosion, step S.-Sediment from erosion has 
higher concentrations of clay, organic matter, nitrogen, 
and phosphorus than the original soil. The enrichment 
coefficient is the ratio of the concentration in the 
sediment to the concentration in the original soil. The 
higliest enrichment coefficients tend to occur with the 
smallest amounts of erosion. Increasing amounts of 
erosion tend to bring the ratio closer to unity [55]. 
This suggests that the sediment from nil erosion has the 
same composition as the soil and that the enrichment is 
the contribution of interrill erosion which is pre- 
dominantly silt and clay. Enriclunent coefficients of two 
are quite common for clay, organic matter, etc., from 
erosion plots. 

The erosion model supplies estimates of the clay 
enrichment coefficients, EC, the pounds of interrill and 
rill erosion, GI and GR, for each compartment in the 
watershed, and the clay fraction, CF, and silt fraction, 
SF, of the surface soil in each zone. The pounds of clay 
in the interrill erosion, PC, is the total clay in the 
sediment minus the clay from rill erosion: 

PC - (CF • EC) (GI + GR)h:F*GR [55] 

and the pounds of sill, PS, is the remainder. 

There are numerous references, see Bailey and 
White (7.)^ to the fact that the different size fractions of 
the soil adsorb chemicals to different degrees. The 
amount of chemical per unit of soil at the surface conies 
from the mean concentration during runoff and equa- 
tion [48]. The fraction of this that is adsorbed to the 
clay is PF * CF/(PF*CF+SF) where PF is the preference 
factor of the chemical for the clay relative to the silt 
when there is none adsorbed on the sand. Therefore, the 
pounds of chemical per pound of clay, SC, is: 

CF PF*CF + SF 

and the remainder of the chemical is associated with the 
sUt. Similar calculations with an additional preference 
factor can be used when the chemical is associated witli 
the sand fraction. 

The total amount of chemical removed in interrill 
erosion is the sum of the product of concentration and 
amount for the clay and silt. 

AI = SC'PC + SS‘PS [57] 

Rill erosion, step 6.-The basic logic used for the loss of 
::hemicals in rill erosion is that all the chemical within 
he rill volume is lost. The area of the rills, EA, as a 
action of the treated area and the average depth of the 


rills, ED, IS provided by the eiosion model. Rills 
probably wider at the soil surface, but for siinplicit; 
rectangular rill shape is assumed. The method 
estimating the width and depth of the rills is ciiide 
this time and doesn’t merit a more complex sha 
consideration. The amount of chemical within llic 
volume is the product of the amount of chemical in t 
compartment, the area of the rills, and the fraction 
the vertical distribution intercepted by the rill. T1 
fraction, Yl, is the integral of a unit distribution fre 
the surface to the bottom of the rilF 

ED 

YI = J* f(x)dx J5i 

0 


When the chemical lias been applied to the surfai 
and leached D inches into the soil, then f (X) is a be 
shaped curve about D and the integral is 



where the minus sign is used when the argument 
negative. A band width of two standard deviation 

2* vH-DC^D, is used to determine how much of th 
chemical distribution the rill depth intercepts. When tli 
chemical is applied by mixing AP inches into the soil an 
then leached D inches, a plug with sharp bouiidariej 
f(x) “ constant, between D and D + AP is assumed. Thi 
creates three limits: 

a. When ED is less than D then 

Yl ^ 0 [60 

b. When ED is greater than D but less than D + A) 
then 

Yl ==: ED/(D + AP) [61] 

c. When ED is greater than D + AP then 

Yl - 1 [62] 

Cultivation that reaches the chemical is assumed to 

create a plug distribution and destroy the previous rills. 
Without cultivation the rill erosion is assumed to occur 
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in the old rills » and the effect of the combined nil depth 
must be used in setting the integral Yl. Summarizing, 
(he amount of chemical lost in rill erosion is 

AR = A*EA-Y1 [63] 

Transport to the channel, step 7 —Runoff water and 
sediment from one xone can cascade over other zones. 
The transformation of watershed into a series of tubes 
was used to describe this feature. The amounts of mnoff 
water, sediments, and any chemical associated with them 
arc accumulated for all the compartments in a tube. 
Currently no allowance is made for loss of the chemical 
over untreated ground by infiltration or adsorption on 
tire soil surface. One estimate of the amount lost by 
infiltration is the product of the infiltration during 
mnoff and the concentration equilibrated with the water 
and sediment of runoff in that compartment. In most 
cases these effects are considered to be small. Since the 
soils could be different in their adsorption and pref- 
erence characteristics as well as texture, the amounts of 
clay and silt arc calculated. The absorption and pref- 
erence factors are weighted for the amount of sediment 
from eacli zone, and the chemical distribution between 
the sediment and water is calculated for the outflow of 
eacli compartment. 

Deposition, step 8. -In some compartments the flow of 
mnoff water will be insufficient to maintain all the 
sediment in suspension -and the excess sediment will be 
deposited, Tlic hydrodynamic characteristics of the 
sediment cause the larger, heavier particles to settle out 
first. In most cases these heavier particles are sand grains 
containing a negligible amount of chemical, It is recog- 
nized that sometimes water-stable aggregates and 
organic-coated grains may cause a significant fraction of 
the chemical load to be deposited with the sand fraction. 
The chemical is distributed according to the weighted 
adsorption coefficient between water and sediment 
before tlie deposition. The amount of deposition is 
satisfied by depositing the sand, silt, and clay in that 
order. The weighted preference factor is used to cal- 
culate the deposition of the chemical with the different 
size fractions. 

Lateral outflow, step 9.— The hydrology model provides 
estimates of the vertical and lateral flow rates for each 
layer in each zone. The compartments at the end of 
tubes are next to the channel, and chemicals in these 
compartments can be carried into the channel by lateral 
flow. Water is also contributed to the outflow of each 
Lube. 

The amount of lateral waterflow is the ratio of the 
lateral flow rate to the total flow rate times the amount 


of drainage The amount of drainage from each layer is 
the amount of infiltration to each layer minus the 
amount required to reach field capacity. Since the 
chemical could be distributed between tsvo layers, the 
fraction in each layer must be calculated. For a bell- 
shaped distribution the fraction above the boundary is 
given by equation [58] where ED is replaced by the 
distance to the boundary. The fraction for a plug 
distribution is the difference between the boundary and 
the top of the plug divided by tlie thickness of the plug. 
The chemical outflow from each layer is then the 
amount of chemical in solution times the outflow 
divided by the amount of water in the soil. 

Mixing for watershed outflow, step lO.-The outflow 
from a watershed is not usually a Well-mixed system. 
Pairs of samples taken at the same time can differ in 
concentration from 9 to 33 percent (77). Some of this 
variability can be attributed to the cascading effect 
symbolized by the tube transformation in this model. 
Thus one estimate of the peak solution concentration in 
the outflow is from the average solution concentrations 
in each tube outflow. The combination of all the tubes 
gives the watershed outflow with a mean solution 
concentration calculated by distributing ail the chemical 
between all the sediment and all the water. 

Location of the remaining chemical, step 11, -The 
presence of a lateral component of waterflow results in 
a diagonal movement of the chemical through the soil 
profile. The total distance moved has been calculated 
from the infiltration. The vertical and horizontal com- 
ponents are calculated from each increment of move- 
ment using the ratio of lateral to vertical flow rates 
weighted for the fraction of the chemical in each layer. 
The peak solution concentration is calculated from 
equations [52] or [54], The deposited chemical is 
added to the chemical already in that compartment or 
creates a new treated area if there is no previous 
chemical, A plug distribution is assumed with a thickness 
of 1 inch plus the depth of the chemical layer already 
there. 

The logic for cultivation effects.-If the. depth of 
cultivation is deeper than half the chemical, it is assumed 
that a plug distribution is created from the surface to the 
depth of tlie cultivation. 

Iteration.-The preceding calculations are performed for 
each compartment in the watershed although some steps 
can be omitted if no chemical is present, After the 
calculations are completed for the entire watershed, the 
results arc totaled and the status of the system is printed 
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out The program then returns to the beginning foi 
another storm event. 

The Nitrate Option 

Nitrate is derived in the soil not only by fertilizer 
application, but also by the mineralization process of 
micro-organisms. Nitrogen in organic matter is converted 
to ammonium, then to nitrite, and finally nitrate. Nitrite 
seldom occurs in significant concentiations and while 
ammonium occurs at significant concentrations more 
frequently, it is moderately adsorbed to the soil. Thus 
the nitrogen problem is usually restricted to nitrate 
which is nonadsorbed and leached easily. Plant uptake, 
denitrification, and immobilization are other mechanisms 
for nitrate removal. 

Since erosion can remove adsorbed organic forms of 
nitrogen and ammonium, but does not significantly 
affect the nitrate movement, the erosion model was 
bypassed and the storm file, from the hydrology model 
normally used by the erosion model, is used directly. 
The elimination of the erosion model also eliminates the 
need to consider compartments within the hydrologic 
zones. Normal mineralization means that the nitrate 
movement in all zones has to be calculated, whether the 
zones received fertilizer or not. The starting date for the 
simulation is the date of the first application of 
fertilizer. The model provides the option for a split 
application of fertilizer to be applied on a second date. 
Table 7 presents an outline of steps used in calculating 
mtrate movement. 

Mineralization, step 1. -Organic nitrogen is mineralized 
to nitrate according to a first order rate equation (45). 
The pounds per acre of nitrate, AN, from mineralization 
during a period of TI days is; 

AN = PN[1.-EXP(RK • TI)] [64] 

where PN is the pounds/acre of potentially mineralizable 
nitrogen in the top layer of soil, and RK in reciprocal 
days is the mineralization rate coefficient. The rate 
coefficient is sensitive to temperatures (^6) and mois- 
ture (47). The Arrhenius equation relating the rate 
coefficient to TM, the average absolute temperature 
(273 + degrees Celcius) during the period, is 

In(RK) - 15.807 “ 6350./TM . [65] 

The Optimum moisture for mineralization appears to be 
near field capacity, and the rate coefficient- decreases 
linearly as the moisture content decreases: 

WC - SW/FC [66] 


Table Z-Outline oj steps for tracing nitrate movement 

1, Mineralization, temperature, and moisture corrections. 

2, Plant uptake, 

3, Fertilizer and cultivation. 

4 Vertical movement and lateral outflow. 

5. Iteration and band completion. 

where WC is the correction for the soil moisture content, 
SW, averaged over the period of mineralization. When 
the moisture content is above field capacity, the conver- 
sion of ammonium to nitrate is inhibited but the 
conversion of organic nitrogen to ammonium is not. 
Short peiiods of ponding are incorporated in the model 
tliroiigii the second infiltration stage. Any ammonium 
buildup is dissipated quickly to nitrate when aerobic 
conditions return. 

Ponding does accentuate denitrification which is not 
accounted for in this model. Considerable study was 
given to tlris process before it was decided to neglect the 
piocess in this version. The presence of a carbon source 
IS necessary for the conversion of nitrate to gases, and 
this would require additional modelling of soluble 
carbon compounds. The same problem was faced and 
the same decision was reached for the process of 
immobilization, the conversion of nitrate into organic 
nitrogen compounds. These are areas for future research. 

Plant uptake, step 2.-The uptake of nitrate by the crop 
IS the mechanism of greatest nitrate loss during the time 
when plants are actively growing. Unpublislied work 
indicates that when the accumulated mineralization plus 
fertilizer does not greatly exceed the amount of nitrogen 
for optimum plant growth (44) only a small amount of 
nitrate remains in the root zone at plant maturity. 

Laboratory studies (5i) show that nitrogen uptake 
increases with the increasing concentration of nitrate in 
solution. Mass flow is the predominant mechanism for 
moving mtrate through the soil to the plant root (2). 
Thus nitrogen uptake should be closely connected to 
evapotranspiration. Evapo transpiration is also related to 
the amount of adsorbing root surface and reflects the 
growth rate of the plant. There is little growth, ET or 
mtrate uptake under cold or drought conditions (32). 
Based on these considerations the nitrate taken up by 
the plant, UP pounds/acre, during the period between 
storms is: 

UP AU • ET/(SW • WD) [67] 

where AU is the pounds of nitrate per acre available for 
uptake, the initial content plus half the mineralized 
amount; ET is the iriches of evapotranspiration from 
each layer weiglited for the distribution of the nitrate 
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within hiyc!^ iiiul the piocliicl of SW and WD gives the 
inclics of walei in tlie soil 

Fertilization and cultivation, step 3,— Fertilizer added 
before a storm is not considered available for uptake 
because uptake fioni split applications of nitrogen is 
small if theie isn't adequate moisture. The applied 
ferlilizcr is available for leaching duiiiig the storm and 
for uptake after the storm. Cultivation is not considered 
directly in the nitrate model, but its effects on water 
movement arc accounted for in the hydrology model. 

Movement, step 4.- The same logic is used for the 
leaching of .nitrate as is used in the general chemical 
nfodel. Infiltration causes vertical movement and lateral 
outflow which is added to the surface runoff. There is 
no nitrate cairied in the suifacc runoff except as 
contribtited by return How, and no loss by erosion or 
deposition is calculated. The model could be modified to 
handle a surface application of nitiogen fertilizer using a 


TESTS FOR 

Carbofuran Experiment at Coshocton, Ohio 

Complete data sets of chemical, crosional, and hydro- 
logical information arc very scarce. The only complete 
set we have been able to compile is for a 2-yoar 
experiment with the pesticide carbofuran on watershed 
No, 113 at (he North Appalachian Fxperimcntal Wiiter- 
shed, Coshocton, OIrio. The data set was developed from 
the miimiseript reporting the iesulls(<5) and supple- 
mented witli information from the imtiiors and the staff 
at the watershed. 


Input to the Hydrological Submodel 

Tables 8 and 9 give the input format for the 
paranielcrs and data. The I’.5-acre watershed was divided 
into three equal zones as shown in figure 1 7. The upper 
zone had an average length of 80 ft and slope of 3 
percent. The middle zone was 110 ft with a 10-percont 
slope and tlic lower zone was 190 ft with 20-percent 
slope. While there is subsurface movement of water in 
litis watershed, none of it comes to lire surface above the 
measuring weir. Corn was grown on the whole watershed 
die 2 years of the experiment, 1971 and 1972. To 
initialize tire hydrologic conditions, the hydrology silb- 
model was operated for the 2 years preceding the 
experiment when the watershed was in meadow. Tem- 
perature and pan evaporation were entered as weekly 


dissolution rate rather than assuming complete dissolu- 
tion in the initial infiltration 

Iteration and band completion, step S.-Tlie preceding 
calculations are repeated for each zone to provide the 
watershed output. All the calculations are repeated for 
each storm event causing the band of nitrate to move 
downward. Since mineralization is assumed to occur m 
tlie top layer, the mineralization subsequent to each 
storm in effect creates a layer that is thicker by the 
amount of movement. When the mineralization rate is 
slow enough, its contribution to the band of nitrate can 
be bypassed and this defines the location of the top of 
the band. This decision is made when the average 
temperature is 5 C or less. The temperature relation 
indicates that the mineralization rate at 5° C is about 10 
percent of optimum. 

Any infiltration during this cool weather (fall, winter 
and spring), causes the band to move further. No explicit 
provisions have been made for frozen conditions, These 
arc presumably handled by the hydrology model. 


VALIDATION 

averages of the daily values, mean degrees Fahrenheit 
and inches. 

The rainfall input was breakpoint values. Each ob- 
servation contained the date, military time, and the 
inches of precipitation since the previous reading. The 
time increment between observations did not exceed 2 
months when there was no rainfall and 30 minutes 
during storm. 

Input to the Erosion Submodel 

The information about cacli compartment needed by 
the erosion submodel is given in table 10. 

In watershed 1 13 all the overland flow cascades to 
the next zone so there is only one stream tube or path, 
The hydrologic variables required are the runoff volume, 
peak runoff rate, and the rainfall intensity for each zone 
in each storm. The slope length and steepness are the 
same as used in the hydrology. In addition, the soil 
texture, 15 percent clay and 60 percent silt, the 
erodibilily factor, 0.48, (he crop factor for corn chang- 
ing from 0.29 to 0.43 with season, and the management 
factor of 0,5, 0.6, and 0.9 for these zones are used, 

Input to the Chemical Submodel 

Tabic IJ lists the information required by tl 
cliemical submodel for each compartment. Only oik 
layer per compartment is necessary because of the 
limited leaching. 
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Figure 17.-^V^LtcIshed 113 divided into three hydrological zones. 


From the hydrology submodel the values of the 
initial soil moisture, infiltration, runoff, and subsurface 
movement are needed. The erosion submodel supplies 
values of interrill and rill erosion, area and depth of nils, 
deposition, and composition of the sediment. The 
chemical characteristics used are given in table 12. In 
1971 the carbofuran was broadcast on the soil surface 
and disced 3 inches into the soil before tlie corn was 
planted. In 1972 the carbofuran was banded 2 inches 
deep with the corn seed. The adsorption coefficient was 
calculated from the concentration on the sedinient and 
in the solution while the breakdown coefficient was 
reported in the manuscript* The banding method in- 
creased the adsorption coefficient threefold and de- 
creased the breakdown coefficient twofold. The disper- 
sion distribution coefficient and the preference factor 
were chosen arbitrarily as ‘‘reasonable” values after 
consulting tire literature, 

Results 

The predicted values of five variables are compared 
with measured values for five storms in table 13. In 1971 
rain started in the evening following application and 
continued the next morning to give the only major 
storm that season. Five* other storms with less than 0.05 
in. of runoff were measured but not predicted by the 
hydrology model. 

In 1972, 26 days passed from the time of application 
until the first storm runoff was measured. Rain totaling 


0.64 in. fell during this period and was treated as initial 
infiltration. Duiing the 25 days between the second and 
third storms, 2.49 in. of ram fell and was treated as part 
of the initial infiltration. EiglU other storms with less 
than 0,08 in. runoff were measured, but not predicted by 
the hydrology submodel before the November?, 1972, 
stoim. 

Table 10,-Ewsion submodel input data for each 
compartment 

Constant watershed characlciistics 

1. Soil credibility factor 

2. Slope steepness, percent 

3. Length of overland flow, ft 

4. Width of compartment, ft 

5. Texture 

Variable watershed characteristics 

6. Crop factor 

7. Management factor 

Storm characteristics 

8. Rainfall factor from storm intensity 

9. Volume of runoff, in 

10. Peak rate of runoff, in/hr 


Table 11. -Information required for the chemical 
submodel 


A, Chemical 


B. 


1. Date of application 

2. Appheationrate, Ibs/a 

3. Depth of mbong, in 

4. Adsorption coefficient, ratio of concentratioms 
Breakdown coefficient, day"^ 


6* Relative preference for size fractions 
Physical 


7. Dispersion distribution factor, in 

8. Field capacity, % 

9. Porosity, % 

10. Texture of soil 


C Hydrological 

1 1 . Date of storm 

12. Infiltration before runoff, in 

13. Infiltration during runoff, m 

14. Runoff, in 

15. Ratio of lateral to vertical flow 

16. Soil moisture at start of storm 

D, Erosional 


17. Interrill erosion, tons 

18. Rill erosion, tons 

19. Fraction of area in rills 

20. Depth of rills, in 

21. Deposition, tons 

22. Texture of sediment 
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Table 12 -Chaiactenstics of catbofuran experiments 
on Coshocton Watershed 1 13 


Date of application 

May 5. 1971 

May 18, 1972 

Amount applied, Ib/a 

4.83 

2.77 

Application method 

Broadcast 

Banded 

Depth of mixing, in ► ♦ » * 

3 

2 

Adsorption coefficient 

Bjcakdown coeflicicnt, day’^ 

1 4 

4.1 

0 015 

0.0074 

Dispeijjion distribution coefficient, 


.5 

5 

Preferciue factor 

10 

10 


The predicted chemical loss and concentration arc 
high. Only preliminaiy ctTorts have been made so far to 
evaluate the cause of the disciepancies fiom the tiiree 
models. Additional results from other chemicals, storms, 
and wateisheds are needed before a thorough evaluation 
can be made 

A moie detailed analysis of the 1971 storm is given in 
table 14. About 64 percent of the chemical loss is 
removed by runoff water in this storm. A similar analysis 
of the 72 storms shows a decrease in the fraction of the 
chemical loss fiom runoff water as the chemical is 


Table 13. -Measured and predicted water, sediment, and carbofuran losses from 
Coshocton Watershed 113 for 5 storms 


Item 



Storm date 




5'6-71 

6-13-72 

6*15-72 

7-10-72 

11-7-72 

Runoff, in 

Meas 

0,21 

0,24 

0.41 

0.20 

0.22 

Pred 

.21 

.46 

.27 

,001 

,15 

Peak flow, cfs 

Meas 

.28 

.85 

3.54 

4,24 

... 

Pred 

,23 

1.15 

.83 

.01 

27 

Plosion, tons/a 

Meas 

1.73 

2.12 

4.49 

.69 

,38 

Pred 

4,27 

9,53 

4.80 

.80 

3.60 

Carbofuran, lbs 

Meas 

.036 

.021 

,046 

.004 

,001 

Pred 

.071 

.15 

.104 

.04 

.013 

Concentration, ppni 

Meas 

.47 

.19 

,22 

.05 

.02 

Pred 

.78 

.56 

.70 

4.24 

.14 


Table 14.— The movement of water, sediment, and carbofuran in each comjxirtment for the 1971 storm 


Item 


Compartment 


Total 

watershed 


1 

2 

3 


Initial infiltration, In 


1.55 

1.99 


Remaining infiltration, in 

99 

.84 

1.38 

... 

Runoff water, in 


.62 

-.36 

... 

Runoff chemical, lb 


,0454 

-.0268 

.0451 

Interiill soil, tons 

1 

.1 

.5 

— 

Interrill chemical, lb 

0002 

.0002 

.0003 

.0007 

Rill soil, tons 

0 

1.1 

4.9 

— 

Rill chemical, lb 

0 

.0046 

.02 

,0246 

Lcachmg depth, in 


.87 

1.03 


Remaining chemical, Ib/a 

4.7 

4.6 

4.8 

— 
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leached deeper into the soil The negative values for 
runoff water and chemical in compartment three oc^ 
curred because the infiltration in this compartment 
exceeded the rainfall and part of the cascading water 
from the second compartment was used. 


Erosion Studies at Treynor, Iowa 

Since additional complete chemical, erosional, and 
hydrological data sets are not presently available, otjier 
data sets liave been used to evaluate parts of the model. 
The hydrology submodel has been evaluated with a 
number of data sets ( 20 ) and will not be discussed here, 
A set of hydrological and erosional data from watershed 
2 at Treynor, Iowa, was obtained from the staff of the 
North Central Watershed Research Center, Columbia, 
Mo., and used to evaluate the combination of the hy- 
drological and erosional submodels. 


Input to the Hydrological Submodel 

A map of the 82.8-acre watershed is shown in figure 
18. The watershed was divided into three zones, 200 ft; 
166 ft, and 133 ft m average overland flow length. It was 
estimated that twenty percent of the upper zone 
cascaded to the middle zone while 15 percent of that 
zone cascaded onto the lower zone. In 1965 com was 
grown, and in the 10-day period from June 25 to July 4 
six major runoff events occurred. To initialize the 
hydrological conditions the hydrology submodel was 
started the previous year when corn was also grown. 


Input to the Erosional Submodel 

To compute the erosion variables the watershed must 
be transformed. The transformation creates five stream 
tubes or paths and seven compartments as shown in 
figure 19. For this silt loam soil the texture is 15 percent 
clay and 65 percent silt, the erodibility factor is 0,32, 
the crop factor for corn is 0.5, and the management 
factor is 0.8 in the upper zone and 0.7 for the other 
zones. 


Results 

Comparisons of the measured and predicted mnoff 
characteristics are given in table 15, for the six storms. 
In general, the predicted values exceed the measured 


values with the gieatest error associated with the largest 
storm. The erosion-deposition chaiactcristics for these 
six storm events arc given m table 16. The only large 
discrepancy between measuied and picdicted sediment 
yield is for the June 29 storm when the hydrologic 
characteristics were also most in error. The detailed 
erosion results of each compartment arc given in table 
17 for the second storm on June 28, 1965. These results 
illustrate the deposition of sediment in the lower 
compartment of the D stream tube 



Figure 18 -Treynor, Iowa Watershed 2. (82.8 acres) 
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Figure 19.— Treynor Watershed W-2 transformed. 
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Table ]3.-RunoJJ chaiactenstics, Treymr Watershed W-2, 6 storm events- 6-27-65 to 

7-1-65 


Date 

Ramfall 


Runoff 


Measured 

Predicted 

Volume 

Peak 

Volume 

Peak 


Jn 

In 

Cfs 

In 

Cfs 

6-27 

0 53 

0.04 

20 

0 16 

51 

6-28 

.62 

.16 

72 

33 

132 

6-28 

1.20 

.68 

152 

92 

204 

6-29 

2.42 

.64 

157 

1.51 

384 

7-1 

.68 

.13 

51 

.12 

39 

7-1 

.25 

.05 

21 

02 

3 


Table 16.- Erosion-deposition characteristics, Treynor Watershed W-2, 6 storm events- 6-27-65 to 7-1-65 


Date 

Rainfall 

amount 

R 

Measured 

sediment 

yield 

Predicted 

sediment 

yield 

Erosion 



Intcrnll 

RiH 



In 


The 

The 

The 

The 

The 

6-27 

0.53 

8.8 

0.51 

0.98 

0.47 

0.61 


6-28 

62 

8,5 

1 84 

1,49 


1.26 

.14 

6-28 

1.20 

31.8 

5.36 

5.74 

1.36 

4.90 

.52 

6-29 

2.42 

26.4 

3.65 

8.21 

1.01 

7.95 

.76 

7-11 

68 

1.6 

.58 

.33 


.30 

.03 

7-1 

25 

.5 

.15 

.05 





Table 1 7. -Erosion results in each compartment for the 2d storm on June 28, 1965 


Compartment Interrlll Deposition 


Amount Area Depth 


Tons Tons Percent In Tons 

A1 26.8 38.0 17.7 0.065 0 

B1 4.7 6.7 17,7 .065 0 

B2 4.1 63.1 24.6 .528 0 

Cl 39.8 192.9 24.6 .248 0 

D1 11.5 55.6 24.6 .248 0 

D2 8.6 24.8 20.0 .170 43.5 

El 17.4 24.8 20.0 .064 0 
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FUTURE PROSPECTS 


As indicated in the discussion of the submodels, the 
hydrology submodel has been tested separately on 
several watersheds and found to produce volumes and 
peaks corresponding reasonably well to the obseived 
data. The coupling of the hydrology and erosional 
submodels has been tested twice, once with the five 
storms on the Coshocton, Ohio watershed and with the 
six storms on the Treynor, Iowa waterslied. The 
combination of all three submodels has been tested once 
with the five storms on the Coshocton, Ohio watershed. 

An adequate test of this model requires complete sets 
of chemical, erosional, and hydrologic data for several 
diverse watersheds. This information is extremely rare 
because it is only lecently that the chemical and physical 
data have been collected smiuUaneously. The only pub- 
lished data that approach the necessary reqiiiiements arc 
the pesticide work on the Coshocton, Ohio watersheds. 
Pesticide experiments recently completed on the Riesel, 
Texas watersheds and currently m progress on the 
Watkinsville, Georgia watershed may provide some 
additional tests of the model 

The best tests of any model come Irom experiments 
designed for that purpose. Therefore, it is hoped that the 
chemists and engineers at all of our watershed centers 
can use this model as a research tool for designing 
experiments and data collection, The experimental 
results will provide tests for this version of the model, as 
well as information for developing better versions. 

These experiments will require an unpredictable 
period of time to complete because of erratic weather. 
In the meantime, the literature is being examined for 
data that miglit serve to test one or two parts of the 
model. A number of papers have been published on the 
chemicals in runoff and erosion from simulated rainfall 
experiments. Some of these may have sufficient data to 
test the combination of erosion and chemical submodels 
for a single compartment. A number of leaching experi- 
ments miglit be used to test the leaching aspect of the 
hydrology and chemical submodel combinations. A few 
of tlie numerous nitrogen experiments might have 
sufficient data to test some aspects of the hydrology and 
nitrate submodel combination. 

It is anticipated that the testing will result in some 
model concepts being revised or replaced. It is hoped 
that tliis model will serve as a stimulus for conducting 
small experiments to test parts of the model; leaching, 
dispersion, rill development, etc. The development of 
the submodels revealed a number of areas where basic 
research is needed to improve our understanding of 


important processes The following few paragraphs 
highlight some of the most obvious problems. 

The path that the water follows through the soil lakes 
on added importance when the chemical aspect is 
considered. The lateral flow of water below the surface 
is extremely difficult to monitor and desciibe. Interflow 
appears to be the best explanation for the highest 
chemical concentration in runoff occurimg in the second 
and third storms of some expeiiments. Snow distribu- 
tion and melting are poorly understood, paiticulaily in 
the piediction of runoff and leaching of chemicals. 
Considerable information and a number of evapotrans- 
piration models have been produced recently, but the 
root growth and water extraction pattern in the soil is 
still inadequately understood. 

The dynamic ot time-dependent aspects of erosion 
take on added importance when chemical adsorption 
and distribution is considered. The detachment, lians- 
poit, and deposition piocesscs depend upon rates as well 
as volumes of lainfall and lunoff. The changing rates of 
watei and soil movement will be icflectcd in time- 
dependent chemical movement. 

Enrichment ratios indicate the need to chaiactcrize 
the sediment movement beyond amount and time. Even 
particle size, while very important, is not sufficient to 
account for the ciiemical transport properties of ag- 
gregates and organic matter. If the concepts of rill and 
inteinll erosion used in this model are shown to be 
important characteristics of chemical loss by erosion, 
then a great deal of researcli will be needed to 
adequately describe such processes as the formation and 
configuration of rills. 

In the chemical area the washoff of pesticides from 
vegetation and the leacliing of nutrients out of the 
vegetation are some of the most poorly understood 
processes. Management practices, such as banding or 
repeated applications, appear to cause greater effects 
than expected. The immobilization and denitrification 
of nitrate continue to be grey areas needing research. 
The chromatographic theory used in this model is being 
tested in the laboratory and in small plots and if found 
to be adequate, considerable work will be needed to 
measure and predict dispersion distribution coefficients. 

The limited target of sloping Corn Belt watersheds 
must be broadened to include all important agricultural 
lands. It is believed that the submodels are adequate for 
most soils. The soils with large cracks are expected to 
cause the most serious problems. Application of the 
model to flat or level lands, such as in parts of the Corn 
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Belt, the delta, and coastal lands, may require alternate 
submodels. The use of large amounts of chemicals in 
these areas makes them a most important challenge. 
Irrigated lands are also flat and receive considerable 
amounts of chemicals. Only parts of the present hy- 


drology and erosion submodels may be needed to model 
the flat and irrigated lands. Salinity has not been 
included because it was felt that there were several 
modeling attempts already underway and thus more 
could be accomplished in these other areas. 
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